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In this paper, Winters & al. (1995) 's available potential energy (APE) framework is 
extended to the fully compressible Navier-Stokes equations, with the aims of: 1) clarifying 
the nature of the energy conversions taking place in turbulent thermally-stratified fluids; 
2) clarifying the role of surface buoyancy fluxes in Munk & Wunsch (1998) 's constraint 
on the mechanical energy sources of stirring required to maintain diapycnal mixing in 
the oceans. The new framework reveals that the observed turbulent rate of increase in 
the background gravitational potential energy GPE r , commonly thought to occur at the 
expenses of the diffusively dissipated APE, actually occurs at the expenses of internal 
energy, as in the laminar case. The APE dissipated by molecular diffusion, on the other 
hand, is found to be converted into IE, similarly as the viscously dissipated kinetic energy 
KE. Turbulent stirring, therefore, does not introduce a new APE / GPE r mechanical- 
to-mechanical energy conversion, but simply enhances the existing IE/GPE r conversion 
rate, in addition to enhancing the viscous dissipation and entropy production rates. This 
in turn implies that molecular diffusion contributes to the dissipation of the available 
mechanical energy ME = APE + KE, along with viscous dissipation. This result has 
important implications for the interpretation of the concepts of mixing efficiency "/mixing 
and flux Richardson number R f , for which new physically-based definitions are proposed 
and contrasted with previous definitions. 

The new framework allows for a more rigorous and general re-derivation from first 
principles of Munk & Wunsch (1998)'s constraint also valid for a non-Boussinesq ocean: 

~ ^—A-OLw — 1 + (1 ~ Olmixing w 

O^jV S2j J ~ Vv r,f orcing — ^ r,f orcing : 

S ttf "/mixing 

where G(KE) is the work rate done by the mechanical forcing, W r j orc i ng is the rate 
of loss of GPE r due to high-latitude cooling, and £ a nonlinearity parameter such that 
£ = 1 for a linear equation of state (as considered by Munk & Wunsch (1998)), but £ < 1 
otherwise. The most important result is that G(APE), the work rate done by the surface 
buoyancy fluxes, must be numerically as large as W r j orc i ng , and therefore as important 
as the mechanical forcing in stirring and driving the oceans. As a consequence, the overall 
mixing efficiency of the oceans is likely to be larger than the value "/mixing = 0.2 presently 
used, thereby possibly eliminating the apparent shortfall in mechanical stirring energy 
that results from using "/mixing = 0.2 in the above formula. 
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1. Introduction 

1.1. Stirring versus mixing in turbulent stratified fluids 

As is well known, stirring by the velocity field greatly enhances the amount of irreversible 
mixing due to molecular diffusion in turbulent stratified fluid flows, as compared with 
the laminar case. A rigorous proof of this result exists for thermally-driven Boussinesq 
fluids for which boundary conditions are either of no-flux or fixed temperature. In that 
case, it is possible to show that 



i.e., the ratio of the entropy production (in the Boussinesq limit) of the stirred state over 
that of the corresponding purely conductive non-stirred state is always greater than unity, 
where T and T c are the temperature of the stirred and conductive states respectively, the 
proof being originally due to Zeldovich (1937), and re-derived by Balmforth & Young 
(2003). The function $ was introduced by Paparella & Young (2002) as a measure of 
the strength of the circulation driven by surface buoyancy fluxes. However, because $ is 
analogous to an average Cox number (the local turbulent effective diffusivity normalised 
by the background diffusivity) e.g. Osborn & Cox (1972); Gregg (1987), it is also 
representative of the amount of turbulent diapycnal mixing taking place in the fluid. 

Reversible stirring and irreversible mixing, e.g., see Eckart (1948), occur in relation 
with physically distinct types of forces at work in the fluid. Stirring works against buoy- 
ancy forces by lifting and pulling relatively heavier and lighter parcels respectively, thus 
causing a reversible conversion between kinetic energy (KE) and available potential 
energy (APE). Mixing, on the other hand, is the byproduct of the work done by the 
generalised thermodynamic forces associated with molecular viscous and diffusive pro- 
cesses that relax the system toward thermodynamic equilibrium, e.g. de Groot & Mazur 
(1962); Kondepudi & Prigogine (1998); Ottinger (2005). Thus, stirring enhances the 
work rate done by the viscous stress against the velocity field, resulting in enhanced 
dissipation of KE into internal energy (IE). Similarly, stirring also enhances the ther- 
mal entropy production rate associated with the heat transfer imposed by the second 
law of thermodynamics, which results in a diathermal effective diffusive heat flux that 
is increased by the ratio (A tur b u ient / Ai am i nar ) 2 (another measure of the Cox number), 
where A tur buient and Ai am i nar refer to the "turbulent" and "laminar" area of a given 
isothermal surface, see Winters & d'Asaro (1996) and Nakamura (1996). In the lam- 
inar regime, the generalised thermodynamic forces associated with molecular diffusion 
are known to cause the conversion of IE into background gravitational potential energy 
(GPE r ). From a thermodynamic viewpoint, it would be natural to expect the stirring to 
enhance the IE / GPE r conversion, but in fact, the existing literature usually accounts 
for the observed turbulent increase in GPE r as the result of a "new" energy conver- 
sion irreversibly converting APE into GPE r . Clarifying this controversial issue is a key 
objective of this paper. 



The most rigorous existing theoretical framework to understand the interactions between 
the different forces at work in a turbulent stratified fluid is probably the available poten- 
tial energy framework introduced by Winters & al. (1995), for being the only one so far 
that rigourously separates reversible effects due to stirring from the irreversible effects due 
to mixing (see also Tseng & Ferziger (2001)). As originally proposed by Lorenz (1955), 
such a framework separates the potential energy PE (i.e., the sum of the gravitational 
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1.2. The modern approach to the energetics of turbulent mixing 
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potential energy GPE and internal energy IE) into its available (APE = AGPE + AIE) 
and non-available (PE r = GPE r + IE r ) components, with the IE component being ne- 
glected for a Boussincsq fluid, the case considered by Winters & al. (f 995). The usefulness 
of such a decomposition stems from the fact that the background reference state is by 
construction only affected by diabatic and/or irreversible processes, so that understand- 
ing how the reference state evolves provides insight into how much mixing takes place in 
the fluid. 

In the case of a freely-decaying turbulent Boussinesq stratified fluid with an equation 
of state linear in temperature, referred to as the L-Boussincsq model thereafter, Winters 
& al. (1995) show that the evolution equations for KE, APE = AGPE, and GPE r 
take the form: 

d(KE) 



dt 

d(APE) 
dt 

d(GPE, 



= —C(KE, APE) - D(KE), (1.2) 
= C(KE, APE) - D(APE), (1.3) 



where C(APE, KE) = -C(KE, APE) is the so-called "buoyancy flux" measuring the re- 
versible conversion between KE and APE, D(APE) is the diffusive dissipation of APE, 
which is related to the dissipation of temperature variance x, e.g. Holloway (1986); 
Zilitinkcvich & al. (2008), while W^mizmg is the rate of change in GPE r induced by 
molecular diffusion, which is commonly decomposed into a laminar W r ^ aminar and tur- 
bulent W r turbulent contribution. All these terms are explicitly defined in Appendix A 
for the L-Boussinesq model, as well for a Boussinesq fluid whose thermal expansion 
increases with temperature, called the NL-Boussinesq model. Appendix B further gen- 
eralises the corresponding expressions for the fully compressible Navier-Stokes equations 
(CNSE thereafter) with an arbitrary nonlinear equation of state (depending on pressure 
and temperature only though). 

Of particular interest in turbulent mixing studies is the behaviour of W r j U rbuient ~ 
the turbulent rate of increase in GPE r — which so far has been mostly discussed in the 
context of the L-Boussinesq model, for which an important result is: 

W rMrhulent = D(APE), (1.5) 

which states the equality between the APE dissipation rate and W r ^ U rbuient - This result 
is important, because from the known properties of D(APE), it makes it clear that 
enhanced diapycnal mixing rates fundamentally require: 1) finite values of APE, since 
D(APE) = when APE = 0; 2) an APE cascade transferring the spectral energy of the 
temperature (density) field to the small scales at which molecular diffusion is the most 
efficient at smoothing out temperature gradients. The discussion of the APE cascade, 
which is closely related to that of the temperature variance, has an extensive literature 
related to explaining the k~ 3 spectra in the so-called buoyancy subrange, both in the 
atmosphere, e.g. Lindborg (2006) and in the oceans, e.g. Holloway (1986); Bouruet- 
Aubertot & al. (1996). Note that because APE is a globally defined scalar quantity, 
speaking of APE cascades requires the introduction of the so-called APE density, noted 
< E> a (x, t) here, for which a spectral description is possible, e.g. Holliday & Mclntyre (1981); 
Roullet & Klein (2009); Molemaker & McWilliams (2009). 

Eqs. (1.2-1.4) exhibit only one type of reversible conversion, namely the "buoyancy 
flux" associated with the APE IKE conversion, and three irreversible conversions, viz., 
D(KE), D(APE), and W r ^ m i X i ng , the first one caused by molecular viscous processes, 
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and the latter two caused by molecular diffusive processes. A primary goal of turbulence 
theory is to understand how the reversible C(APE, KE) conversion and irreversible 
D(KE), D(APE), Wr^mixmg are all inter-related. In this paper, the focus will be on tur- 
bulent diffusive mixing, for the understanding of viscous dissipation constitutes somehow 
a separate issue with its own problems, e.g. Gregg (1987). The nature of these links is 
usually explored by estimating the energy budget of a turbulent mixing event, defined 
here as a period of intense mixing preceded and followed by laminar conditions, for which 
there is a huge literature of observational, theoretical, and numerical studies. Integrating 
the above energy equations over the duration of the turbulent mixing event thus yields: 



AKE = -C{KE, APE) - D(KE), (1.6) 



AAPE = C{KE, APE) - D(APE), (1.7) 

AGPE r W r,mixing — r, turbulent ~l~ r, laminar i (^-^) 

where A(.) and the overbar denote respectively the net variation and the time-integral 
of a quantity over the mixing event. Summing the KE and APE equation yields the 
important "available" mechanical energy equation: 



AKE + AAPE = ~[D(KE) + D(APE)} < 0, (1.9) 

which states that the total "available" mechanical energy ME = KE + APE undergoes a 
net decrease over the mixing event as the result of the viscous and diffusive dissipation of 
KE and APE respectively. A schematic of the APE dissipation process, which provides 
a diffusive route to KE dissipation, is illustrated in Fig. 1. 

1.3. Measures of mixing efficiency in turbulent stratified fluids 

Eq. (1.9) makes it clear that turbulent diapycnal mixing (through D(APE)) participates 
in the total dissipation of available mechanical energy ME = KE+APE. Since D(APE) 
is non-zero only if APE is non-zero, turbulent diapycnal mixing therefore requires having 
as much of ME in the form of APE as possible. The classical concept of "mixing effi- 
ciency" , reviewed below, seeks to provide a number quantifying the ability of a particular 
turbulent mixing event in dissipating ME — KE + APE preferentially diffusively rather 
than viscously. From a theoretical viewpoint, it is useful to separate turbulent mixing 
events into two main archetypal categories, corresponding to the two cases where ME 
is initially entirely either in KE or APE form. These two cases are treated separately, 
before providing a synthesis addressing the general case. 

At a fundamental level, quantifying the mixing efficiency of a turbulent mixing event 
requires two numbers, one to measure how much of ME is viscously dissipated, the other 
to measure how much of ME is dissipated by turbulent mixing. While everybody seems 
to agree that D(KE) is the natural measure of viscous dissipation, it is the buoyancy 
flux C(APE, KE), rather than D(APE), that has been historically thought to be the 
relevant measure of how much of ME is dissipated by turbulent mixing, since it is 
the term in Eq. (1.6) that seems to be removing KE along with viscous dissipation. For 
mechanically-driven turbulent mixing events, defined here as being such that AAPE = 
and AME — AKE, the efficiency of mixing has been classically quantified by means of 
two important numbers. The first one is the so-called flux Richardson number Rf, defined 
by Linden (1979) as: "the fraction of the change in available kinetic energy which appears 
as the potential energy of the stratification" , mathematically defined as: 



= C(KE,APE) = C(KE,APE) 
f " \AKE\ C{KE,APE) + D{KEY 
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e.g., Osborn (1980), and the second one is the so-called "mixing efficiency": 



Rf C(KE,APE) 
7™ s - 1 _ Rf ~ ■ U-llj 

It is now recognised, however, that the buoyancy flux represents only an indirect measure 
of irreversible mixing, since it physically represents a reversible conversion between KE 
and APE, while furthermore appearing to be difficult to interpret empirically, e.g. see 
Barry & al. (2001) and references therein. Recognising this difficulty, Caulfield & Peltier 
(2000) and Staquet (2000) effectively suggested to replace C(KE, APE) by a more direct 
measure of irreversible mixing in the above definitions of Rf and ^mixing- Since turbulent 
diapycnal mixing is often diagnosed empirically from measuring the net changes in GPE r 
over a mixing event, e.g. McE wan (1983a,b); Barry & al. (2001); Dalziel & al (2008), 
a natural choice is to use W ' r turbulent as a direct measure of irreversible mixing, which 
leads to: 

j^GPEr- _ W r ,turbulent ^ ^ 

S W r ,turbulent + D(KE) 

GPE r ll f _ VV r.turbulent -. „s 

^mixing ~ j _ R GPE r ~ D ( KE ) ' <■ ' > 

From a theoretical viewpoint, these definitions are justified from the fact that in the 
L-Boussinesq model, the following equalities hold: 



C{APE, KE) = D(APE) = W rMrbulent , (1.14) 

as follows from Eqs. (1.6) and (1.7), combined with Eq. (1.5), when AAPE = 0. The 
modified flux Richardson number Rj PEr coincides — for a suitably defined time interval 
— with the cumulative mixing efficiency £ c introduced by Caulfield and Peltier (2000), 
as well as with the generalised flux Richardson number R^ defined by Staquet (2000), in 

GPE 

which our 7 mia;i ^g is also denoted by 7&. 

Although Eqs. (1.12) and (1.13) are consistent with the traditional buoyancy-flux-based 
definitions of Rf and ^mixing in the context of the L-Boussincsq model, such definitions 
overlook the fact that Eq. (1.14) is not valid in the more general context of the fully 
compressible Navier- Stokes equations, for which the ratio 

>. Wr, turbulent /-. _\ 

^= D(APE) (L15) 

is in general less than one, and even sometimes negative, for water or seawater. For this 
reason, it appears that Rf and jmixing should in fact be defined in terms of D(APE), 

not W r turbulent, viz., 



Rf— = D{APE) (1.16) 
1 D(APE) + D(KE) 



DAPE D(APE) 

Imixing DIKE) ' \ ') 

which we call the dissipation flux Richardson number, and dissipation mixing efficiency 
respectively, to distinguish them from their predecessors. In our opinion, Rf lAPE and 
Imixing as defined by Eqs. (1.16) and (1.17) are really the ones that are truly consistent 
with the properties assumed to be attached to those numbers. Most notably, Eq. (1.16) 
is the only way to define a flux Richardson number that is guaranteed to lie within 
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the interval [0,1], since neither C(KE, APE) nor W r ,turbulent can be ascertained to be 
positive under all circumstances. Since Eqs. (1.12) and (1.13) are still likely to be used 
in the future owing to their practical interest, it is useful to provide conversion rules 
between the GPE r and D(APE)-based definitions of Rf and ^mixing, viz., 

t uDAPE 

GPE r t^DAPE r>GPE r 2_f (~\-[fl) 

Imixing S Imixingi lx -f 1 — (1 — £ t )R DAPE ' 

These formula require knowledge of the nonlincarity parameter £, which measures the 
importance of nonlinear effects associated with the equation of state, see Tailleux (2009) 
for details about this. The often-cited canonical value for mechanically-driven turbulent 
mixing is j m ixing ~ 0.2, which appears to date back from Osborn (1980), e.g. Peltier & 
Caulfield (2003). 

The second case of interest, namely buoyancy-driven turbulent mixing, is defined here 
as being such that AKE = and AME = AAPE, as occurs in relation with the 
so-called Rayleigh- Taylor instability for instance. Eqs. (1.6) and (1.7) lead to: 



C(KE,APE) = -D(KE) < (1.19) 



D(APE) = C{KE, APE) - AAPE = \ AAPE\ - \C(KE, APE)\. (1.20) 

Eq. (1.19) reveals that the buoyancy flux is negative this time, and that it represents 
the fraction of ME that is lost to viscous dissipation, not diffusive dissipation. This 
establishes, if needed, that the buoyancy flux should not be systematically interpreted 
as a measure of irreversible diffusive mixing. Since Linden (1979)'s above definition for 
the flux Richardson number does not really make sense for Rayleigh- Taylor instability, 
an alternative definition is called for. The most natural definition, in our opinion, is as 
the fraction of ME dissipated by irreversible diffusive mixing, viz., 



_ -AAPE + C{KE, APE) _ \C{KE,APE)\ 
f ~ -AAPE ~ 1 \AAPE\ ' 

which according to Eqs. (1.6) and (1.7), is equivalent to: 



(1.21) 



_ DiAPE, 
3 D(APE) + D(KE) 



with the corresponding value of Jmixing'- 



7™ = -5*- = (1.23) 
1 9 1 - R f D{KE) ' 

which are identical to Eqs. (1.16) and (1.17). The above results make it possible, therefore, 
to use R EAPE and J^i^ing as definitions for the flux Richardson number and mixing 
efficiency that make sense for all possible types of turbulent mixing events. 

At this point, a note about terminology seems to be warranted, since in the case 
of the Rayleigh- Taylor instability, it is Rf that is referred to as the mixing efficiency 
by some authors, e.g. Linden and Redondo (1991); Dalziel & al (2008), rather than 
^mixing- Physically, this seems more logical, since Rf is always comprised within the 
interval [0, 1], whereas ^mixing is not. Interestingly, Oakey (1982) appears to be the first 
to define ^mixing as a: "mixing coefficient representing the ratio of potential energy to 
kinetic energy dissipation" . For this reason, it would seem more appropriate and less 
ambiguous to refer to ^mixing as the "dissipations ratio" . Unfortunately, it is not always 
clear in the literature which quantity the widely used term "mixing efficiency" refers to, 
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as it has been used so far to refer to both Rf and "/mixing- In order to avoid ambiguities, 

DAPE 

/mixing 5 



the remaining part of the paper only make use of the quantities Rj JAPE and j^ AJ t 



which for simplicity are denoted Rf and "/mixing respectively. 

As a side note, it seems important to point out that Raylcigh- Taylor instability has the 
peculiar property that AGPE r , max , the maximum possible increase in GPE r achieved 
for the fully homogenised state, is only half the initial amount of APE, e.g. Linden and 
Redondo (1991); Dalziel & al (2008) (at least when £ = 1, i.e., in the context of the 
L-Boussinesq model). Physically, it means that less than 50 % of the initial APE can 
actually contribute to turbulent diapycnal mixing, and hence that at least 50 % of it 
must be eventually viscously dissipated. As a result, one has the following constraints: 



_ D(APE) _ V, turbulent 1 , , 

Hf - \AAPE\ ~ \AAPE\ 2 [ ' 

"/mixing ^5 ^ ~£J2 ^ ^ (^'^^) 

Experimentally, Linden and Redondo (1991) reported values of Rf « 0.3 ("/mixing = 
3/7 w 0.43), while Dalziel & al (2008) reported experiments in which the maximum 
possible value Rf « 0.5 ("/mixing ~ 1) was reached. Owing to the peculiarity of the 
Rayleigh- Taylor instability, however, one should refrain from concluding that "/mixing = 1 
or Rf = 0.5 represent the maximum possible values for "/mixing and Rf in turbulent 
stratified fluids. To reach definite and general conclusions about "/mixing and Rf, more 
general examples of buoyancy-driven turbulent mixing events should be studied. It would 
be of interest, for instance, to study the mixing efficiency of a modified Rayleigh- Taylor 
instability such that the unstable stratification occupies only half or less of the spatial 
domain, so that AGPE rmax \AAPE\. In this case, all of the initial APE could in 
principle be dissipated by molecular diffusion, which would correspond to the limits 
Rf = 1 and "/mixing = +oo. Of course, such limits cannot be reached, as it is impossible 
to prevent part of the APE to be converted into KE, part of which will necessarily 
be dissipated viscously, but they are nevertheless important in suggesting that values of 
"/mixing > 1 can in principle be reached, which sets an interesting goal for future research. 

1.4. On the nature of D (APE) and W r ,turbuient 
Of fundamental importance to understand the physics of turbulent diapycnal mixing are 
the nature and type of the energy conversions associated with D(APE) and W T . turbulent- 
So far, it seems fair to say that these two energy conversions have been regarded as 
essentially being one and the same, based on the exact equality W r j U rbulent — D(APE) 
occurring in the L-Boussinesq model, suggesting that molecular diffusion irreversibly con- 
verts APE into GPE r , e.g. Winters & al. (1995). Such an interpretation appears to be 
now widely accepted, e.g. Caulfield & Peltier (2000), Peltier & Caulfield (2003), Munk 
& Wunsch (1998), Huang (2004), Thorpe (2005) among many others. The main char- 
acteristic of this view, schematically illustrated in the top panel of Fig. 2, is to disregard 
the possibility that the turbulent increase of GPE r might be due the enhancement of 
the IE /GPE r conversion rate by the stirring. In other words, the current view assumes 
that the work involved in the turbulent increase of GPE r is done by the stirring against 
buoyancy forces, not by the generalised thermodynamic forces responsible for entropy 
production and the IE/GPE r conversion. At the same time, the current view seems to 
accept that stirring enhances entropy production. But from classical thermodynamics, 
this seems possible only if the work rate done by the generalised thermodynamic forces 
is also enhanced, which in turn should imply an enhanced IE/GPE r conversion. 
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In order to determine whether the turbulent increase of GPE r could be accounted 
for by a stirring-enhanced IE/GPE r conversion rate, rather than by the irreversible 
conversion of APE into GPE r , it seems useful to point out that the validity of Winters 
& al. (1995)'s interpretation seems to rely crucially on D(APE) and W r j U rbulent being 
exactly identical, not only mathematically (as is the case in the L-Boussinesq model) 
but also physically. Here, two quantities are defined as being physically equal if they 
remain mathematically equal in more accurate models of fluid flows — closer to physical 
"truth" in some sense — such as CNSE for instance. Indeed, only a physical equality can 
define a physically valid energy conversion, as we hope the reader will agree. As shown 
in Appendix B, however, which extends Winters & al. (1995) results to the CNSE, 
the equality D(APE) = W r ,turbuient is found to be a serendipitous feature of the L- 
Boussinesq model, which at best is only a good approximation, the general result being 
that the ratio 

£ Wr .turbulent 

4 - D(APE) ■ [L2b) 

usually lies within the interval — oo < £ < 1 for water or seawater, and that it strongly 
depends on the nonlinear character of the equation of state. Whether there exists fluids 
allowing for £ > 1 is not known yet. An important result is that it appears to be perfectly 
possible for GPE r to decrease as the result of turbulent mixing, in contrast to what is 
often stated in the literature. This case, which corresponds to £ < 0, was in fact previously 
identified and discussed by the late Nick Fofonoff in a series of little known papers, see 
Fofonoff (1962, 1998, 2001). For this reason, the case £ < shall be subsequently 
referred to as the Fofonoff regime, while the more commonly studied case for which 
W r ^turbuient > shall be referred to as the classical regime. 

The lack of physical equality between D(APE) and W r ^ U rbulent makes Winters & 
al. (1995) 's interpretation very unlikely, and gives strong credence to the idea that 
W r . turbulent actually correspond to a stirring-enhanced IE/GPE r conversion rate. If 
so, what about D(APE)? In order to shed light on the issue of APE dissipation, it is 
useful to recall some well known properties of thermodynamic transformations associated 
with the following problem: Assuming that the potential energy PE = GPE + IE of 
a stratified fluid increases by AE, how is AE split between AGPE and AIE? Here, 
standard thermodynamics tells us that the answer depends on whether AE is added 
reversibly or irreversibly to PE. Thus, if AE is added reversibly to PE (i.e., without 
entropy change, and for a nearly incompressible fluid), then: 

AGPE , AIE , , . 

-AE- A£ « 1 ( L2? ) 

while if AE is added irreversibly (i.e., with an increase in entropy), then: 

AGPE _ AIE , nn . 

-AE- <<1 > ~AE~ ~ (L28) 

i.e., the opposite. These results, therefore, suggest that when molecular diffusion converts 
APE into PE r , the dissipated APE must nearly entirely go into IE r , not GPE r , in 
contrast to what is usually assumed (The demonstration of Eqs. (1.27) and (1-28) is 
omitted for brevity, but this follows from the results of Appendix B.) It follows that 
what the equality D(APE) — W r j U rbulent of the L-Boussinesq actually states is the 
equality of the APE /IE and IE/GPE r conversion rates (or more generally, for real 
fluids, the correlation between the two rates), not that D(APE) and W r ^ U rbuient are of 
the same type. Physically, the two conversion rates W r j U rbulent and D(APE) appear to 
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Figure 1. Idealised depiction of the diffusive route for kinetic energy dissipation. (I) represents 
the laminar state possessing initially no AGPE and AIE, but some amount of KE. (II) represents 
the state obtained by the reversible adiabatic conversion of some kinetic energy into APE, which 
increases APE but leaves the background GPE r and IE r unchanged; (III) represents the state 
obtained by letting the horizontal part of molecular diffusion smooth out the isothermal surfaces 
until all the APE in (II) has been converted into background PE r = GPE r + IE r . 



be fundamentally correlated because they are both controlled by molecular diffusion and 
the spectral distribution of APE, as will be made clear later in the text. 



1.5. Internal Energy or Internal Energies? 

In the new interpretation proposed above, internal energy is destroyed by the IE/GPE r 
conversion at the turbulent rate W r ^urbuient, while being created by the APE conversion 
at the turbulent rate D(APE). Could it be possible, therefore, for the dissipated APE to 
be eventually converted into GPE r , not by the direct APE/GPE r conversion route pro- 
posed by Winters & al. (1995), as this was ruled out by thermodynamic considerations, 
but indirectly by transiting through the IE reservoir? 

As shown in Appendix B, the answer to the above question is found to be negative, 
because it turns out that the kind of IE which APE is dissipated into appears to be dif- 
ferent from the kind of IE being converted into GPE r . Specifically, Appendix B shows 
that IE is indeed best regarded as the sum of distinct sub-reservoirs. In this paper, 
three such sub-reservoirs are introduced: the available internal energy {AIE) , the exergy 
(IE exergy ), and the dead internal energy (IEq). Physically, this decomposition parallels 
the following temperature decomposition: T{x,y, z,t) — T' \x,y, z,t) + T r (z,t) + To(t), 
where T (t) is the equivalent thermodynamic equilibrium temperature of the system, 
T r {z,t) is Lorenz's reference vertical temperature profile, and T" (x,y, z,t) the residual. 
Physically, AIE is the internal energy component of Lorenz (1955)'s APE, while IEq 
and IE exergy are the internal energy associated with the equivalent thermodynamic equi- 
librium temperature T and vertical temperature stratification T, r respectively. The idea 
behind this decomposition can be traced back to Gibbs (1878), the concept of exergy 
being common in the thermodynamic engineering literature, e.g. Bejan, A. (1997). See 
also Marquet (1991) for an application of exergy in the context of atmospheric avail- 
able energetics. A full review of existing ideas related to the present ones is beyond the 
scope of this paper, as the engineering literature about available energetics and exergy 
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A) Standard Interpretation of Eq. (1.5) 
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B) New Interpretation of Eq. (1.5) 
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Figure 2. (A) Predicted energy changes for an hypothetical turbulent mixing event under the 
assumption that the diffusively dissipated APE is irreversibly converted into GPE r ; (B) Same 
as in A under the assumption that the diffusively dissipated APE is irreversibly converted into 
IEo, as the viscously dissipated KE. In both cases, the net energy changes in KE, GPE r , APE, 
and IE are the same. The only predicted differences concern the subcomponents of the internal 
energy IEo and IE exergy . 



is considerable. The way it works is encompassed in the following equations: 

d{KE) 



dt 

d(APE) 



dt 

d(GPE r ) 



= —C(KE,APE) — D{KE), (1.29) 
= C(KE, APE) - D(APE), (1.30) 



dt 

d(IE Q ) 



dt 

d(I Eexergy ) 



^^r,mixing ^^r,laminar turbulent-) 

D(KE) + D(APE) = D total , (1.32) 



^^r,mixing — ^^r,laminar ^^r, turbulent • (1.33) 



dt 

In this model, the first three equations are just a rewriting of Eqs. (1.2) and (1.4), so that 
the main novelty is associated with Eq. (1.32) and (1.33). Physically, Eq. (1.32) states 
that the viscous and diffusive dissipation processes D(KE) and D(APE) mostly affect 
To but not T r , while Eq. (1.33) states that the IE/GPE r conversion reduces IE exergy as 
well as smoothes out T r . The empirical verification of the validity of the above equations 
is the main topic of Section 2. 
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1.6. Link with the ocean heat engine controversy 

In the oceans, turbulent diapycnal mixing is a crucial process, as it is required to trans- 
port heat from the surface equatorial regions down to the depths cooled by high-latitude 
deep water formation. In the traditional picture found in most oceanography textbooks, 
turbulent diapycnal mixing and deep water formation are usually described as part of 
the buoyancy-driven component of the large-scale ocean circulation responsible for the 
oceanic poleward transport of heat, often called the meridional overturning circulation 
(MOC). Physically, the MOC is often equated with the longitudinally-averaged circu- 
lation taking place in the latitudinal/vertical plane. The possible dependence of the 
buoyancy-driven circulation on mechanical forcing, which one might expect in a system 
as nonlinear as the oceans, has been usually ignored. However, the idea of a buoyancy- 
driven circulation unaffected by mechanical forcing physically makes sense only if one 
can establish that the mechanical stirring required to sustain turbulent diapycnal mixing 
is driven by surface buoyancy fluxes. Munk & Wunsch (1998) (MW98 thereafter) ques- 
tioned this view, and argued instead that turbulent diapycnal mixing must in fact be 
primarily driven by the wind and tides, and hence that the buoyancy-driven circulation 
must in fact be mechanically-controlled. Moreover, MW98 analysed the GPE budget of 
the oceans to derive the following constraint: 

G(KE) = ^ r ' /orcmg , (1.34) 

^mixing 

linking the work rate G(KE) done by the mechanical sources of stirring, the rate W r j- orc i n(l 
at which high-latitude cooling depletes GPE r , and the oceanic bulk "mixing efficiency" 
^mixing (or more accurately, the dissipations ratio, as argued previously). Physically, Eq. 
(1.34) states that the fraction ^mixing of G(KE) has to be expanded into turbulent mix- 
ing to raise GPE r at the same rate W r j orc i ng at which it is lost. By using the values 
W r j orcing « 0.4 TW and -/mixing = 0.2^ MW98 concluded that G(KE) = 0(2 TW) is 
approximately required to sustain the observed oceanic rates of turbulent diapycnal mix- 
ing. This result caused much stirring in the ocean community, because the wind supplies 
only about 1 TW, leaving an apparent shortfall of 1 TW to close the energy budget. This 
led MW98 to argue that the only plausible candidate to account for the missing stirring 
should be the tides, spawning a considerable research effort over the past 10 years on the 
issue of tidal mixing. 

Although MW98's arguments have been echoed favourably within the ocean commu- 
nity, e.g. Huang (2004); Paparellafe Young (2002); Wunsch & Ferrari (2004); Kuhlbrodt 
(2007); Nycander & al. (2007), it remains unclear why the surface buoyancy fluxes should 
not be important in stirring and driving the oceans, given that the work rate done by 
the surface buoyancy fluxes, as measured by the APE production rate, was previously 
estimated by Oort & al. (1994) to be G(APE) = 1.2 ± 0.7 TW and hence comparable 
in importance to the work rate done by the mechanical forcing. In their paper, MW98 
rather summarily dismissed Oort & al. (1994) 's results by contending that the so-called 
Sandstrom (1908)'s "theorem" requires that G(APE) be negligible, and hence that the 
buoyancy forcing cannot produce any significant work in the oceans, but given the highly 
controversial nature of Sandstrom (1908) 's paper in the physical oceanography commu- 
nity, and its apparent refutation by Jeffreys (1925), it would seem important to have a 
more solid physical basis to make any definitive statements about G(APE). Note that 
Sandstrom's paper was recently translated by Kuhlbrodt (2008), who argues that Sand- 
strom's did not initially formulate his results as a theorem, but rather as an inference. 
Ascertaining whether G(APE) is large or small is obviously crucial to determine whether 
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the MOC is effectively driven by the turbulent mixing powered by the winds and tides, as 
argued by MW98, or whether it is in fact predominantly buoyancy-driven, as it appears 
to be possible if G(APE) is as large as predicted by Oort & al. (1994). Further clari- 
fication is also needed to understand the possible importance of some effects neglected 
by MW98, such as those due to a nonlinear equation of state, which Gnanadesikan & al. 
(2005) argue lead to a significant underestimation of G(KE), or due to entrainmcnt ef- 
fects, which Hughes & Griffiths (2006) argue lead to a possible significant overestimation 
of G(KE). 

1.7. Purpose and organisation of the paper 

The primary objective of this paper is to clarify the nature of the energy conversions 
taking place in turbulent stratified fluids, with the aim of clarifying the underlying as- 
sumptions entering MW98's energy constraint Eq. (1.34). The backbone of the paper are 
the theoretical derivations presented in the appendices A and B, which provide a rigor- 
ous theoretical support to understand the links between stirring and irreversible mixing 
in mechanically and thermodynamically-forced thermally-stratified fluids. Appendix A 
offers a new derivation of Winters & al. (1995) 's framework, which is further extended 
to the case of a Boussinesq fluid with a thermal expansion increasing with temperature. 
Appendix B is a further extension to the case of a fully compressible thermally-stratified 
fluid, in which the decomposition of internal energy into three distinct sub-reservoirs is 
presented. Section 2 seeks to illustrate the differences between D(APE) and W rt turbuient 
using a number of different viewpoints, and examine some of its consequences, in the 
context of freely decaying turbulence. Section 3 revisits the issues pertaining to MW98's 
energy constraint. Section 4 offers a summary and discussion of the results. 



2. A new view of turbulent mixing energetics in freely decaying 
stratified turbulence 

2.1. Boussinesq versus N on- Boussinesq energetics 

As mentioned above, a central point of this paper is to argue that irreversible energy 
conversions in turbulent stratified fluids are best understood if internal energy is not 
regarded as a single energy reservoir, but as the sum of at least three distinct sub- 
reservoirs. Obviously, these nuances are lost in the traditional Boussinesq description of 
turbulent fluids, since the latter lacks an explicit representation of internal energy, let 
alone of its three sub-reservoirs. This does not mean that the Boussinesq approximation 
is necessarily inaccurate or incomplete, but rather that the definitive interpretation of 
its energy conversions requires to be checked against the understanding gained from the 
study of the fully compressible Navier-Stokes equations. Such a study was carried out, 
with the results reported in Appendix B. In our approach, successive refinements of the 
energy conversions were sought, starting from the KE / APE / PE r system for which the 
number of energy conversions is limited and unambiguous. The second step was to split 
PE r into its GPE r and IE r components; the third step was to further split IE r into its 
exergy IE r — IEq and dead IEq components. Finally, the last step was to split APE into 
its AIE and AGPE components. These successive refinements are illustrated in Fig. 10 
in Appendix B. An important outcome of the analysis is that the structure and form 
of the KE I APE jGPE r equations (Eqs. (1.2-1.4)) obtained for the L-Boussinesq model 
turn out to be more generally valid for a fully compressible thermally-stratified fluid, so 
that one still has: 

^3 = -C(KE,APE) - D(KE), (2.1) 
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(II) HORIZONTAL MIXING 

Figure 3. Idealised depiction of the numerical experimental protocol used to construct Fig. 4, 
as well as underlying the method for constructing Figs. 5 and 6. In the top panel, a piece of 
stratified fluid is cut into pieces of equal mass that are numbered from 1 to Ntot, where N to t is 
the total number of parcels. A random permutation is generated as a way to shuffle the parcels 
randomly and adiabatically, in order to mimic the stirring process. In the bottom panel, all the 
parcels lying at the same level are homogenised to the same temperature by conserving the total 
energy of the system, which mimics the horizontal mixing step illustrated in Fig. 1. 



dt 

d(GPE r ) 



d(APE) = C ^ KE ^ ApE ^ _ D ( ApE ^ ( 2 2) 



W^r,mixing — W^r,turbulent ~t~ Wr, laminar ■ (^-^) 



^ ''remixing 

It can be shown, however, that the explicit expressions for C(KE,APE), D(KE), and 
D(APE) differ between the two sets of equations, see Appendices A and B for the 
details of these differences. Based on the numerical simulations detailed in the following, 
the most important point is probably that D(APE) appears to be relatively unaffected 
by the details of the equation of state, in contrast to W r ^ m i X i ng , which suggests that 
the L-Boussinesq model is able to accurately represent the irreversible diffusive mixing 
associated with D(APE). Moreover, since the internal energy contribution to APE is 
usually small for a nearly incompressible fluid, it also follows that the L-Boussinesq model 
should also be able to capture the time-averaged properties of C(KE, APE), since the 
latter is the difference of two terms expected to be accurately represented by the L- 
Boussinesq model based on the APE equation. The L-Boussinesq model, however, will in 
general fail to correctly capture the behaviour of GPE r , unless the approximation of a 
linear equation of state is accurate enough, as seems to be the case for compositionally- 
stratificd flows for instance, e.g. Dalziel & al (2008). The above properties help to 
rationalise why the L-Boussinesq model appears to perform as well as it often does. 

Being re-assured that there are no fundamental structural differences between the 
energetics of the KE/APE/GPE r system in the Boussinesq and compressible NSE, the 
next step is to clarify the link with internal energy. One of the main results of this paper, 
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(a) (b) 




Diffusively dissipated SE ~ A ('E r -IE ) 

Figure 4. (a) The increase of AGPE as a function of the stirring energy SE (see text for details). 
Each point represents a different stratification shuffled by a different random permutation. The 
continuous line represents the straight-line of equation AAGPE — SE which would describe the 
energetics of the stirring process if AIE could be neglected; (b) the change of GPE r as a function 
of the stirring energy SE dissipated by diffusive mixing; the dotted line is the straight line of 
equation AGPE r = Diffusive dissipated SE which would describe the energetics of turbulent 
mixing if the irreversible conversion AGPE — > GPE r existed; (c) The change in dead internal 
energy IEo as a function of the diffusively dissipated stirring energy SE. The dashed line is 
the straight-line of equation AIEq = diffusively dissipated SE; the figure shows a near perfect 
correlation; (d) The change in GPE r as a function of the exergy change. The dashed line is 
the straight line of equation AGPE r = —AIEexergy- The figure shows again a near perfect 
correlation. 



derived in Appendix B, are the following evolution equations for the dead and exergy 
components of internal energy: 

d ( IE °) _ D{KE) + D(APE), (2.4) 



dt 

d(IE exergy ) 



t, ■- ''r,rmxingi 

which were obtained by neglecting terms scaling as 0(aP/ (pC p )), for some values of a, P, 
p, and C p typical of the domain considered, where a is the thermal expansion coefficient, 
P is the pressure, p is the density, and C p is the heat capacity at constant pressure. 
The important point is that such a parameter is very small for nearly incompressible 
fluids. For seawater, for instance, typical values encountered in laboratory experiments 
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done at atmospheric pressure are a = 2.10~ 4 K~ 1 , P = 10 5 P a , C p = 4.10 3 J.K _1 .kg _1 , 
p = 10 3 m 3 .kg _1 , which yield aP/(pC p ) = 5.10 -6 . In the deep oceans, this value can 
increase up to O(10~ 3 ), but this is still very small. Eqs. (2.4) and (2.5) confirm that 
D(APE) and D(KE) are fundamentally similar dissipative processes, in that they both 
convert APE and KE into dead internal energy, while also confirming that W r ^ m i X i ng 
represent a conversion between IE exergy and GPE r both in the laminar and turbulent 
cases. 

2.2. Analysis of idealised turbulent mixing events 

To gain insight into the differences between D(APE) and W r turbulent ■> the energy bud- 
get of an hypothetical turbulent mixing event associated with shear flow instability is 
examined in the light of Eqs. (2.1-2.5). Typically, such events can be assumed to evolve 
from laminar conditions with no APE. Once the instability is triggered, APE starts 
to increase and oscillate until the instability subsides and the fluid re-laminarises, at 
which point APE returns to zero. The mixing event causes the shear flow to lose a cer- 
tain amount of kinetic energy \AKE\, as well as GPE r to increase by a certain amount 
AGPE r , as the result of the partial smoothing out of the mean vertical temperature gra- 
dient by molecular diffusion. Integrating Eqs. (2.1-2.5) over the duration of the mixing 
event yields: 



AKE = -C(KE, APE) - D(KE), (2.6) 



AAPE = = C{KE, APE) - D(APE), (2.7) 

^GP E r — W r,turbulent ~i~ r,laminar •> (^*^) 



AIEq = D(APE) + D(KE), (2.9) 

AI E exer gy \W r^lamin ar ~t~ I 7 ^ r.turbulent] ; (^'-^) 

where A(.) and the overbar denote a quantity's net change over the time interval and 
its time- integrated value respectively. From an observational viewpoint, energy conver- 
sion terms such as C(KE, APE) are difficult to measure directly; moreover, the results 
can be ambiguous, e.g. Barry & al. (2001) and references therein. As a result, energy 
conversions are probably best inferred from measuring changes in the different energy 
reservoirs, as this appears to be easier to do accurately. Multiple possible inferences arise, 
however, if IE variations are not separated into their IE exergy and IE components. Fig. 
2 illustrates this point for an hypothetical turbulent mixing experiment with hypothetical 
plausible numbers, by showing that a given observed net change in IE of +0.79 units can 
potentially be explained — in absence of any knowledge about the respective variations 
in IE n and IE exergy — as either due to the conversion of 0.8 unit of KE into IE minus 
the conversion of 0.01 unit of IE exergy into GPE r , or by the conversion of 0.8 unit of 
KE into IEq plus the conversion of 0.2 unit of APE into IE® minus the conversion of 
0.21 unit of IE exergy into GPE r . Although the first interpretation is the one implicit in 
Winters & al. (1995) and currently favoured in the literature, it is not possible, based 
on energy conservation alone, to reject the second interpretation. In fact, the only way 
to discriminate between the two interpretations requires separately measuring IEo and 
IE exergy variations, as only then are the two interpretations mutually exclusive. 

2.3. An idealised numerical experimental protocol to test the two interpretations 

To compute I Eg and IE exergy , a knowledge of the temperature field and of the Gibbs 
function for the fluid considered (see Feistel (2003) in the case of water or seawater) 
is in principle sufficient. It is hoped, therefore, that the present study can stimulate 
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laboratory measurements of IEq and IE exerg y , in order to provide experimental support 
(or refutation, as the case may be), for the present ideas. In the meantime, numerical 
methods are probably the only way to assess the accuracy of the two key formula Eqs. 
(2.4) and (2.5), which physically argue that: 1) the diffusively dissipated APE is nearly 
entirely converted into IE ; and 2) that GPE r variations are nearly entirely accounted 
for by corresponding variations in IE exergy . 

To prove our point, energetically consistent idealised mixing events are constructed and 
studied numerically. The procedure is as follows. One starts from a piece of thermally 
stratified fluid initially lying in its Lorenz (1955)'s reference state in a two-dimensional 
container with a flat bottom, vertical walls, and a free surface exposed at constant atmo- 
spheric pressure at its top. The fluid is then discretised on a rectangular array of dimen- 
sion N x x N z into discrete fluid elements having all the same mass Am = pAxAz, where 
x and z are the horizontal and vertical coordinates respectively, as illustrated in Fig. 3. 
The initial stratification has a vertically-dependent temperature profile T(x, P) = T r (P) 
regarded as a function of horizontal position x and pressure P. Thousands of idealised 
mixing simulations are then generated according to the following procedure: 

(a) Initialisation of the reference stratification. The initial stratification is discretised 
as T it k = T(xi,Pk) = T r (P k ), with x t = (i - l)Ax, i = 1,...N X and P k = P min + 
(k — l)gAm, k = 1, . . . N z , where P m in and T r (Pk), k = 1, . . . N z are random generated 
numbers such that T m i n < T r (Pk) ^ T max that have been re-ordered in the vertical to 
create a statically stable stratification, for randomly generated T m i n , T max , and P m in- 

(b) Random stirring of the fluid parcels. The fluid parcels are then numbered from 
1 to N = N x x N z , and randomly shuffled by generating a random perturbation of N 
elements, such that each parcel conserves its entropy in the re- arrangement. Such a step is 
intended to mimic the adiabatic stirring of the parcels associated with the KE — ► APE 
conversion. The random stirring of the fluid parcels requires an external amount of energy 

- called the stirring energy SE — which is diagnosed by computing the difference in 
potential energy between the shuffled state and initial state, i.e., 

SE = (GPE + IE) shuffled - (GPE + IE) initial . (2.11) 

The latter computation requires a knowledge of the thermodynamic properties of the 
fluid parcels. In this paper, such properties were estimated from the Gibbs function for 
seawater of Feistel (2003) by specifying a constant value of salinity. Thermodynamic 
properties such as internal energy, enthalpy, density, entropy, chemical potential, speed 
of sound, thermal expansion, haline contraction, and several others, are easily estimated 
by computing partial derivatives with respect to temperature, pressure, salinity, or any 
combination thereof, of the Gibbs function. The stirring energy SE is none other than 
Lorenz (1955) 's APE of the shuffled state. Since the stirring leaves the background 
potential energy unaffected, the energetics of the random shuffling is given by: 

AAPE = AAGPE + AAIE = SE, (2.12) 

AGPE r = AIE r = 0, (2.13) 

where Eq. (2.12) states that the stirring energy SE is entirely converted into APE, 
while Eq. (2.13) expresses the result that being a purely adiabatic process, the stirring 
leaves the background reference quantities unaltered. Fig. 4 (a) depicts AAGPE as a 
function of SE for thousands of experiments, all appearing as one particular point on 
the plot. According to this figure, AAGPE approximates AAPE within about 10%. 
This illustrates the point that for adiabatic processes, APE is well approximated by its 
gravitational potential energy component, as expressed by Eq. (1.27). 



Energetics and thermodynamics of turbulent molecular diffusive mixing 17 



(c) Isobaric irreversible mixing of the fluid parcels. In the last step, all the fluid parcels 
lying in the same layer are mixed uniformly to the same temperature, by assuming an 
isobaric process that conserves the total enthalpy of each layer. Such a process converts 
a fraction qSE of the APE into the background PE r , according to: 

AAPE = AAGPE + AAIE = -qSE (2.14) 

AG PE r + AIE r = qSE (2.15) 

where < q ^ 1. The factor q is needed here because mixing each layer uniformly does 
not necessarily lead to a statically stable stratification; when this happens, the resulting 
stratification still contains some APE — (1 — q)SE associated with the static instability, 
so that q — 1 only when the mixed density profile is statically stable. 
The change in AGPE r resulting from the irreversible mixing step is depicted as a function 
of the diffusively dissipated stirring energy qSE in the panel (b) of Fig. (4). If the 
stirring energy were entirely dissipated into GPE r , as is classically assumed, then all 
points should lie on the line of equation AGPE r = qSE appearing as the dashed line in 
the figure. Even though such a relation appears to work well in a number of cases, the 
vast majority of the simulated points corresponds to cases where AGPE r is significantly 
smaller than qSE, and even often negative as expected in the Fofonoff regime discussed 
above. On the other hand, if one plots AIE as a function of the diffusively dissipated 
stirring energy qSE, as well as AGPE r as a function of the exergy change AIE exergy = 
—A(IE r — IEo), as done in panels (c) and (d) of Fig. (4) respectively, then a visually 
near perfect correlation in both cases is obtained. This is consistent with the following 
relations: 

AIEq w qSE, (2.16) 

AGPE r « -A(IE r - IE ), (2.17) 

and hence in agreement with the approximate Eqs. (2.4) and (2.5). Eq. (2.16) empirically 
verifies the Eq. (1.28). QED. 

2.4. Numerical estimates of B, W r . m ixing and D( APE) as a function of APE 

Having clarified the nature of the net energy conversions occurring in idealised mixing 
events, the following turns to the estimation of the turbulent rates of the three important 
conversion terms B, W r , mixing, and D(APE), which are the main three terms affected 
by molecular diffusion in the fully compressible Navier-Stokes equations, where B is 
the work of expansion/contraction. As pointed out in the introduction, enhanced rates 
fundamentally arise from turbulent fluids possessing large amounts of small-scale APE. 
For this reason, this paragraph seeks to understand how the values of B, W r , m ixing and 
D(APE) are controlled by the magnitude of APE. 

We first focus on the work of expansion/contraction B, which takes the following form: 

f aP f aP f P DP 

B= / —\7 ■ (k P C p \7T) dV + / — pedV- / —~^dV, (2.18) 

Jv P^P Jv P^P Jv P c s Dt 

obtained by regarding p as a function of temperature and pressure. The part of B affected 
by molecular diffusion is the first term in the right-hand side of Eq. (2.18), and the one 
under focus here. The second and third term in the r.h.s. of Eq. (2.18) are respectively 
caused by the work of expansion due to the viscous dissipation Joule heating, and to the 
adiabatic work of expansion/contraction. The study of these two terms is beyond the 
scope of this paper. 

The diffusive part of B was estimated numerically for thousands of randomly generated 
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Figure 5. (Top panels) The work of expansion/contraction B normalised by its laminar value 
(obtained for APE=0) as a function of a normalised APE for a particular stratification corre- 
sponding to the classical regime, with the right panel being a blow-up of the left panel. (Bottom 
panels) Same as above figure, for the same temperature stratification, but taken at a mean 
pressure of 50 dbar instead of atmospheric pressure, which is sufficient to put the system in the 
Fofonoff regime. The figures show that although B is usually negative in every case, it is never- 
theless positive for small values of APE in the classical regime, as expected from L-Boussinesq 
theory. The normalisation constant APE max corresponds to the overall maximum of APE for 
all experiments. 



stratifications, similarly as in the previous paragraph, using a standard finite difference 
discretisation of the molecular diffusion operator. Unlike in the previous paragraph, how- 
ever, all the randomly generated stratifications were computed from only two different 
reference states pertaining to the classical and Fofonoff regimes respectively, the results 
being depicted in the top and bottom panels of Fig. 5 (with the right panels providing a 
blow-up of the left panels). The main result here is that finite values of APE can make 
the diffusive part of B negative and considerably larger by several orders of magnitude 
than in the laminar APE = case. This result is important, because it is in stark con- 
trast to what is usually assumed for nearly incompressible fluids at low Mach numbers. 
From Fig. 5, it is tempting to conclude that there exists a well-defined relationship be- 
tween the diffusive part of B and APE, but in fact, the curve B = B(APE) is more 
likely to represent the maximum value achievable by B for a given value of APE. Indeed, 
it is important to realize that a given value of APE can correspond to widely different 
spectral distributions of the temperature field. In the present case, it turns out that the 
random generator used tended to generate temperature fields with maximum power at 
small scales, which in turns tend to maximise the value of B for a particular value of 
APE. For the same value of APE, smoother stratifications exist with values of B lying 
in between the x-axis and the empirical curve B = B(APE), the latter being expected 
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Figure 6. (Left panels) The rate of change of GPE r normalised by its laminar value as a function 
of normalised APE, in the classical regime (top panel) as well as in the Fofonoff regime (bottom 
panel). The stratification is identical to that of Fig. 5. (Right panels) The rate of diffusive 
dissipation of APE normalised by W r ,mixing laminar value, as a function of a normalised APE, 
in the classical regime (top panel) , as well as for the Fofonoff regime (bottom panel) . The figure 
illustrates the fact that if the former can be regarded as a good proxy for the latter in the 
classical regime, as is usually assumed, this is clearly not the case in the Fofonoff regime. The 
two figures also illustrate the fact that the former always underestimate the latter for a thermally 
stratified fluid, so that observed values of mixing-efficiencies obtained from measuring GPE r 
variations are necessarily lower-bounds for actual mixing efficiencies. 



to depend on the numerical grid resolution employed. Nevertheless, Fig. 5 raises the in- 
teresting question of whether the empirical curve B = B(APE) could in fact describe 
the behaviour of the fully developed turbulent regime, an issue that could be explored 
using direct numerical simulations of turbulence. 

The two remaining quantities of interest are W rym i X i ng and D(APE), which were nu- 
merically estimated from the following expressions derived in Appendix B: 

W r , mixing = f ^-V ■ (k P C p VT) dV (2.19) 

Jy Pr-Gpr 

D(APE) = j Tr T T X7 ■ (npGpVT) dV = J n P C p X7T ■ (^-jr^j dV. (2.20) 

As for B, these two quantities were evaluated for thousands of randomly generated strat- 
ifications as functions of APE, starting from the same reference states as before. The 
results for W r , mixing are depicted in the left panels of Fig. (6), and the results for D( APE) 
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in the right panels of the same Figure, with the top and bottom panels corresponding 
to the classical and Fofonoff regimes respectively. The purpose of the comparison is to 
demonstrate that whereas there exist stratifications for which the two rates D(APE) 
and W rjm i X i n g are nearly identical (top panel, classical regime), as is expected from the 
classical literature about turbulent stratified mixing, it is also very easy to construct 
specific cases occurring in the oceans for which the two rates become of different signs 
(bottom panel, Fofonoff regime). The other important result is the relative insensitivity of 
D(APE) to the nonlinear character of the equation of state compared to W r ^ mix i ng , sug- 
gesting that the use of the L-Boussinesq model can still accurately describe the KE/APE 
interactions even for strongly nonlinear equations of states, although it would fail to do a 
good job of simulating the evolution of GPE r outside the linear equation of state regime. 
This also suggests that the L-Boussinesq should be adequate enough to study the mix- 
ing efficiency of turbulent mixing events over a wide range of circumstances, provided 
that by mixing efficiency one means the quantity ^mixing = D(APE) / D(KE), and not 
Jmixing = Wr^urbuient I ' D(K E) . Finally, we also experimentally verified (not shown) that 
D(APE) is well approximated by the quantity W r ^mixing — B, as is expected when AIE 
is only a small fraction of APE. 

2.5. Synthesis 

The energetics of freely decaying turbulence is summarised in Fig. 7 for the classical 
(top panel) and Fofonoff (middle panel) regimes, with the bottom panel attempting 
further synthesis by combining AIE and AGPE into a single reservoir for APE, and 
the two regimes into a single diagram. Doing so makes the bottom panel of Fig. (7) 
basically identical to the Boussinesq energy flowchart depicted in the bottom panel of 
Fig. 2. Interestingly, the middle panel suggests that the Fofonoff regime may differ from 
the extensively studied classical regime in several fundamental ways. Indeed, whereas 
both W and B act as net sinks of KE in the classical regime, it appears possible in 
the Fofonoff regime for some fraction of the KE dissipated into AIE to be recycled 
back to KE. This is reminiscent of the positive feedback on the turbulent kinetic energy 
discussed by Fofonoff (1998, 2001), who suggested that such a feedback would enhance 
turbulent mixing and hence speed up the return to the classical regime after sufficient 
reduction of the vertical temperature gradient. If real, such a mechanism would be very 
important to study and understand, as potentially providing a limiting process on the 
maximum value achievable by the buoyancy frequency, with important implications for 
numerical ocean models parameterisations. In his papers, however, Fofonoff envisioned 
the positive feedback on turbulent KE as being associated with the conversion of GPE r 
into AGPE, but this goes against the findings of this paper arguing that GPE r can 
only be exchanged with the exergy reservoir. Fofonoff's feedback mechanism was also 
criticised by McDougall & al. (2003) on different grounds. While the present results 
do not necessarily rule out Fofonoff's feedback mechanism, they suggest that the latter 
probably does not work as originally envisioned by Fofonoff, if it works at all (McDougall 
& al. (2003) 's arguments are not really conclusive either, as they implicitly rely on the 
existence of the APE/GPE r conversion). In any case, the issue seems to deserve more 
attention, given that many places in the oceans appear to fall into Fofonoff's regime. 

3. Forced/dissipated balances in the oceans 

3.1. A new approach to the mechanical energy balance in the oceans 

Prior to re- visiting MW98's energy constraint Eq. (1.34), we start by establishing a 
number of important results for mechanically- and thermodynamically-forced thermally- 
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Figure 7. The energetics of freely decaying turbulence for the classical regime (top panel), the 
Fofonoff regime (middle) panel, and a synthesis of both regimes obtained by subsuming AGPE 
and AIE into APE alone. Note the similarity of the energetics in the lower panel and that of 
the re-interpreted Boussinesq energetics of the lower panels of Fig. 2. 



stratified fluids, based on the results derived in Appendices A and B. The main modifi- 
cations brought about by the mechanical and thermodynamical forcing is the apparition 
of forcing terms, i.e., terms involving the external forcing, in the evolution equations for 
KE, APE, and GPE r as follows: 

d(KE) = _ C ^ KE ^ ApE ^ _ D ( KE ^ + G{KE), (3.1) 



dt 

d(APE) 
dt 

d{GPE r ) 
dt 



= C(KE, APE) - D(APE) + G(APE), (3.2) 

^ T r .turbulent ~t~ ^^r.laminar forcing: ■ (^-^) 



^Vr, mixing 

where G(KE) is the work rate done by the external mechanical forcing, G(APE) is 
the work rate done by the buoyancy forcing, and W r j orc i ng is rate of change of GPE r 
(usually a loss, hence the assumed sign convention) due to the buoyancy forcing. The 
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Figure 8. Energy flowchart for a mechanically- and buoyancy-driven thermally-stratified fluid, 
where Q ne t = Qheating — Qcooiing- At leading order, the "dynamics" (the KE/APE/IEq 
system) is decoupled from the "thermodynamics" (the IE exergy /GPE r system). The dynam- 
ics/thermodynamic coupling occurs through the correlation between D(APE) and W r ,mixing, 
as well as through the correlation between G(APE) and W r , forcing- 



resulting energy transfers are illustrated in the energy flowchart depicted in Fig. 8. This 
figure shows that at leading order, the "Dynamics" — associated with the reservoirs 
KE/APE/IEq — is decoupled from the "Thermodynamics" — associated with the 
GPE r / IE exergy energy reservoirs. Indirect coupling occurs, however, from the fact that 
D(APE) and W r , turbulent on the one hand, and G(APE) and W r j orc i ng on the other 
hand, are strongly correlated to each other. 

3.1.1. Link between G(APE) and W r j orc i ng 

Unlike in the L-Boussincsq model, G(APE) and W r j orc i ng differ from each other in a 
real compressible fluid, for the same reasons that D(APE) differs from W r _turbuient , as 
is apparent from their exact formula given by Eqs. (B 23) and (B 38) in Appendix B: 

G(APE) = [ T Z Tr K pC p \7T ■ ndS, (3.4) 
Js T 



Vjorcing = ~ f "' ^ ^ Kp C p VT ■ ndS. (3.5) 
J S Pr^pr 

In order to understand by how much W r j orc i ng differs from G(APE) in a real fluid, it 
is useful to expand T as a Taylor series around P = P r , i.e., 

T = T(P a ) =T r + T r (P a -P r ) + -.- 

where T r — a r T r /(p r C pr ) is the adiabatic lapse rate, e.g. Feistel (2003). As a result 

T - T r _ T r (P r - P a ) | _ a r (P r - P a ) 

Inserting Eq.(3.6) into Eq. (3.4) reveals that G(APE) and W r j orc i ng are in fact equal 
at leading order. In that case, therefore, the equality between G(APE) and W r j orc i ng 
that exactly holds in the L-Boussinesq model appears to be a much better approximation 
than the corresponding equality between D(APE) and W r ,turbuient- For this reason, we 
shall neglect the differences between G(APE) and W r j or cing in the following. 
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3.1.2. Steady-state Mechanical energy balance 

Under steady-state conditions, summing Eqs. (3.1) and (3.2) yields: 

G{KE) + G(APE) = D(APE) + D(KE), (3.7) 

which simply states that in a steady-state, the production of mechanical energy by the 
wind and buoyancy forcing is balanced by the viscous and diffusive dissipations of KE 
and APE respectively. Fig. 9 schematically illustrates how the wind and buoyancy forcing 
can both contribute to the creation of APE. The main novelty here is to make it clear that 
D(APE) is a "true" dissipation mechanism, i.e., one that degrades mechanical energy 
into internal energy (as does viscous dissipation), not one converting mechanical energy 
into another form of mechanical energy (i.e., GPE r ). This suggests regarding D(APE) + 
D(KE) as the total dissipation of available mechanical energy ME = KE + APE. 

In contrast, most studies of ocean energetics of the past decade have tended to subsume 
the APE production and dissipation terms into the single term B = G(APE) — D(APE), 
in which case Eq. (3.7) becomes: 

G(KE) + B = D(KE). (3.8) 

The problem in writing the mechanical energy balance under this form is that it erro- 
neously suggests that B, rather than G(APE), is the work rate done by surface buoyancy 
fluxes, and that viscous dissipation is the only form of mechanical energy dissipation. For 
instance, Wang & Huang (2005) estimated B = (9(1.5 GW) in the oceans, in the context 
of the L-Boussinesq model, which is about three orders of magnitude less than the work 
rate done by the wind and tides. In one of the most recent review about ocean energetics 
by Kuhlbrodt (2007), it is Wang & Huang (2005) 's estimate for B that is presented as 
the work rate done by surface buoyancy fluxes, while Oort & al. (1994)'s previous result 
for G(APE) is regrettably omitted. 

In the L-Boussinesq model, B takes the particular form B = Kg((T) t0 p — (T) bottom), 
where (T) top and (T) bottom are the area- integrated surface and bottom temperature re- 
spectively, so that in the absence of mechanical forcing, Eq. (3.8) becomes: 

Kg {(T) top - (T) bottom ) = D(KE). (3.9) 

In a study addressing the issue of horizontal convection, recently reviewed by Hughes 
& Griffiths (2008), Paparella & Young (2002) proved that the left-hand side of (3.9) 
must be bounded by n times some finite constant when the fluid is forced by a sur- 
face temperature condition, with no-normal flux applying everywhere else. This result is 
now commonly referred to as the "anti-turbulence theorem" , for the bound implies that 
D(KE) must vanish in the "inviscid" limit (used here to mean both vanishing molecular 
viscosity and diffusivity) , thus violating the so-called "zeroth law of turbulence" , an em- 
pirical law grounded in many observations showing that the viscous dissipation of KE in 
homogeneous turbulent fluid flows remain finite and independent of molecular viscosity 
as the Reynolds number is increased indefinitely. 

3.1.3. Actual implications of the anti-turbulence theorem 

As shown by Wang & Huang (2005), Paparella & Young (2002) 's bound suggests 
that the oceanic viscous dissipation D(KE) would be less than 1.5 GW in absence of 
mechanical forcing. Since this value is several orders of magnitude than observed oceanic 
values of D{KE), the result demonstrates that mechanical forcing is essential to account 
for the latter. By itself, however, the result says nothing about whether mechanical 
forcing is also essential to account for the observed turbulent rates of diapycnal mixing, 
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R a 10 b 2.10 6 3.10 6 4.10 6 5.10 6 
$ 6.2 6.9 7.5 7.9 8.3 

Table 1 . Values of $ as a function of the Rayleigh number R a reproduced from Fig. 4 of 
Paparella & Young (2002) . Values of R a appropriate to the oceans are of the order 

R a = O(10 20 ). 



since, as far we are aware, the values of Rf and ^mixing for horizontal convection have 
never been determined before. Indeed, the anti-turbulence theorem only imposes that 
the difference B = G(APE) - D(APE) be small, but this docs not forbid G(APE) and 
D(APE) to be individually very large. In fact, the current theoretical and numerical 
evidence suggest that G(APE) and D(APE) increase with the Rayleigh number R a . 
Indeed, this is suggested by Paparella & Young (2002) 's numerical experiments, which 
show the function $ (given by Eq. 1.1), which we interpreted as a measure of the cox 
number 0(K t /k) (and hence of D(APE)) to increase with R a = gaATH 3 /(vk) as 
tabulated in Table 3.1.3. Although such values appear to be much smaller than observed 
O(10 2 — 10 3 ) Cox numbers, they also correspond to Rayleigh numbers that are about 
13 to 14 orders of magnitude smaller than occuring in the oceans, leaving open the 
possibility for $ to be possibly much larger, possibly as large as encountered in the 
oceans. In a related study, Siggers & al. (2004) derived a bound on <I> (which they 
related to a horizontal Nussclt number), which does not exclude the possibility that 
horizontal convection, on its own, could support a north-south heat transport of the 
observed magnitude. A further discussion of the physics of horizontal convection based 
on laboratory experiments is provided by Mullarney & al. (2004). In summary, while the 
anti-turbulence theorem demonstrates the need for mechanical forcing to account for the 
observed values of kinetic energy dissipation, the question of whether mechanical forcing 
is needed to sustain diapycnal mixing rates and a north-south heat transport of the 
observed magnitude is still largely open. In particular, it is important to point out that 
although the anti-turbulence theorem rules out the possibility of elevated values of kinetic 
energy dissipation in absence of mechanical forcing, it does not rule out the possibility 
of elevated values of diapycnal mixing rates. In this respect, Paparella & Young (2002)'s 
suggestion that horizontal convection should be regarded as "non-turbulent" appears 
somewhat misleading. 

3.1.4. B acts- of-the- envelope estimate of G(APE) for the world oceans 

In order to make progress on the above issue, it is essential to determine how large 
G(APE) can be in the oceans. As mentioned above, Oort & al. (1994) inferred G(APE) = 
1.2 ± 0.7 TW from observations, and hence to be nearly as large as the work rate done 
by the wind, but this estimate was questioned by MW98 on the basis of Sandstrom's 
"theorem" . A possible source of error in Oort & al. (1994)'s is its reliance on the so-called 
Lorenz approximation, which is often said to overestimate G(APE), e.g. Huang (1998). 

In fact, the simplest method to convince oneself that G(APE) must be large in the 
oceans comes from the result that G(APE) « W r j orc i ng established previously, which 
states that if the rate of decrease of GPE r due to the buoyancy forcing is large, so must it 
be the case for G(APE). Since MW98 inferred W r j orc i ng w 0.4 TW, one can immediately 
conclude G(APE) w 0.4 TW, which it turns out is close to the lower bound of Oort 
& al. (1994) 's estimate, consistent with the idea that the method used by the latter 
should overestimate G(APE). This immediately establishes that MW98's assumption 
that G(APE) is small is inconsistent with their assumption that W r j orc i ng is large. It 
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Figure 9. Idealised depictions of mechanically-driven (left panel) and buoyancy-driven (right 
panel) creation of APE. I) A wind blowing at the surface of a two-layer fluid causes the tilt 
of the layer interface, resulting in a net C(KE, APE) > conversion. II) Localised cooling at 
high-latitudes sets the density of a fraction of the upper layer to that of the bottom layer, also 
inducing a tilt in the layer interface. The return of the interface to equilibrium conditions (flat 
interface) results in a net C(APE,KE) > conversion. 




also establishes that Sandstrom's "theorem", whatever it means, cannot say anything 
meaningful about G(APE). 

An independent way to estimate G(APE) is by using the exact formula for G(APE) 
or Wrjorcmg given by Eqs. (3.4) and (3.5) recalled above. In these formulae, T r is the 
temperature that a surface parcel would have if lifted adiabatically to its reference level. 
Since the oceans are on average heated and cooled where they are the warmest and 
coolest respectively, the parcels' reference level will be on average close to the surface 
in the warm regions, but much deeper in the cold regions. Eq. (3.4) must therefore be 
dominated by surface cooling. Using the near equality between G(APE) and W r j orc i ng , 
we take as our estimate for G(APE) the expression: 

G(APE) « (^p-) QcooUng (3.10) 

\Pr<spr / coo u n g 

with a value of (a r P r / (p r Cp r )) C ooiing appropriate to the regions of cooling. Using the val- 
ues a = HAQ-^K- 1 , P r = 2000 dbar = 2.10 7 Pa, p r = 10 3 kg, C pr = 4.10 3 J.kg-^K" 1 , 
and QcooUng = 2 PW yields: 

G(APE) = 5 ' 1 "~ 5 * 2 -}f x 2.10 15 W - 0.5 TW 
10 3 x 4.10 3 

which is very close to Munk & Wunsch (1998)'s estimate for W r j or cing, an d is consistent 
with Oort & al. (1994)'s lower bound for G(APE). The large value of G(APE) suggests 
that buoyancy forcing can actively participate in maintaining turbulent diapycnal mixing 
in the oceans, which should be reflected by a higher value of ^mixing than the value 
^mixing — 0.2 currently used in the literature about the subject. 

3.2. A new look at the GPE r balance and Munk & Wunsch (1998) 's theory 

Having clarified the "available mechanical energy balance" , we now turn to the GPE r 
budget, with the aim of elucidating the assumptions underlying MW98's constraint on 
the energy requirement for sustaining diapycnal mixing in the oceans. In a steady-state, 
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the GPE r budget given by Eq. (3.3) becomes: 

mixing — ^ r ^turbulent ~t~ ^rjamifiar — ^^r, forcing ; (^'H) 

where the explicit expressions for W r ^ m ixing and W r j orc i ng are given at leading order by 
the following expressions: 



W^r,mixing 



[ K pC p \7T-v(^-)dV, (3.12) 



W r ,mixing 



W rJorcing ^ f (^-) K P C p \7T-ndS, (3.13) 

is \Pr<^pr J 

which are valid for a fully compressible thermally-stratified ocean. Note that for the 
L-Boussinesq model, implicitly considered by MW98, W r ^ m ixing can be rewritten as: 

/ K\\Vz r \\ 2 a r ^-dV w / KrpoN^dV (3.14) 

by using the definition of turbulent diapycnal diffusivity of Winters & al. (1995) for K~t- 
As Eq. (3.14) is exactly the expression used by MW98, this establishes that MW98's 
analysis actually pertains to the GPE r budget, not the GPE budget, and that their 
results should logically follow from Eq. (3.11). 

In order to show that this is indeed the case, simply use the definition "/mixing = 
D(AP E) / D(K E) in combination with the mechanical energy balance to express D(APE) 
in terms of the total mechanical energy input G(APE) + G(KE) as follows: 

D(APE) = 7 """" g [G(APE) + G(KE)] = R f [G{APE) + G(KE)], (3.15) 

1 "T "/mixing 

where Rf — "/mixing I '(1 + Imixing) is the dissipation flux Richardson number defined 
in the introduction. Now, writing W rttU rbuient — t]D{APE) as proposed in this paper to 
account for a nonlinear equation of state, neglecting W r ^ am i nar compared to W r j U rbuienti 
and using the result that W r j orc i ng « G(APE) demonstrated previously, Eq. (3.11) 
becomes: 

^D(APE) w G(APE). (3.16) 
The desired result is obtained by combining Eqs. (3.15) and (3.16), which yields: 

G{RE) ^ l + (l-t)lm l x l ng G{ApE) = G (APE) . (3.17) 

<H mixing K^f 

This formula generalises MW98's Eq. (1.34) to account for a nonlinear equation of state, 
the effects of which being contained in the single parameter £ < 1. It is easily seen that 
MW98's Eq. (1.34) is recovered simply by setting £ = 1 in Eq. (3.17), using the result 
that G(APE) w W r j orc i ng . This formula is further extended and discussed in Tailleux & 
Rouleau (2009), in the context of idealized experiments of mechanically-stirred horizontal 
convection. 

The confirmation that G(APE) should actually be nearly as large as the work rate done 
by the mechanical forcing makes it possible for buoyancy forcing to drive possibly a very 
large fraction of the oceanic turbulent diapycnal mixing, which should be reflected in an 
appropriate value for "/mixing- As noted in the introduction, buoyancy- driven turbulent 
mixing is usually significantly more efficient than mechanically-driven turbulent mixing, 
suggesting that a value of "/mixing significantly larger than the value "/mixing — 0.2 should 
be used in Eq. (3.17). Likewise, the fact that the nonlinear equation of state for seawater is 
found to strongly affect changes in GPE r , as showed in the previous section, suggests that 
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^mixing — 0.2 mixing — 0.5 ^mixing — 1 

£ = 1 2TW 0.8TW 0.4 TW 

£ = 0.5 4.4TW 2 TW 1.2TW 

Table 2. Mechanical energy requirements on G(KE) depending on different assumed values for 
the dissipations ratio ^mixing and nonlinearity parameter £, as computed from Eq. (3.17) using 
G(APE) = 0.4 TW, in line with Munk & Wunsch ( 1998) 's assumptions. 



a value of £ between and 1 should be used (note that £ cannot be negative if a steady- 
state is to exist). A detailed discussion of which values should actually be used for "/mixing 
and £ is beyond the scope of this paper, however, as much more needs to be understood 
about buoyancy-driven and mechanically-driven turbulent mixing in a non-Boussinesq 
stratified fluid before one may become confident enough to speculate on the "right" 
values. In order to fix ideas, however, it is useful to compute the energy requirement 
on turbulent mixing predicted by Eq. (3.17) for plausible values of £ and -) m ixing, as 
reported in Table 3.2, under the assumption that W r j orc i ng « G(APE) w 0.4 TW. 

As expected, decreasing £ at fixed ^mixing increases the requirement on G(KE). This 
is consistent with Gnanadesikan & al. (2005) 's conclusions that cabelling (i.e., the con- 
traction upon mixing stemming from the nonlinear character of the equation of state for 
seawater) increases the requirement on G(KE). Likewise, increasing ^mixing at fixed £ 
decreases the requirement on G{KE). This is consistent with Hughes & Griffiths (2006) 's 
argument that the requirement on G(KE) can be decreased if the entrainment of am- 
bient water by the sinking cold plumes is accounted for. Indeed, taking into account 
entrainment effects is equivalent to increasing j m i X ing, since entrainment is physically 
associated with buoyancy-driven turbulent mixing, as far as we understand the issue. 
Note that decreasing £ seems to require increasing jmixing if one accepts the idea that 
no more than 2TW is available from mechanical energy sources to stir the oceans. Al- 
ternatively, one could also perhaps question the assumption that the oceans are truly in 
a steady state. 

4. Summary and conclusions 

In this paper, we extended Winters & al. (1995) APE framework to the fully com- 
pressible Navier-Stokes equations, with the aims of: 1) clarifying the nature of the energy 
conversion taking place in turbulent stratified fluids; 2) clarifying the role of the surface 
buoyancy fluxes in Munk & Wunsch (1998) 's constraint on the mechanical sources of 
stirring required to sustain diapycnal mixing in the oceans. The most important results 
are that the well-known turbulent increase in background GPE ri commonly thought to 
occur at the expense of the diffusively dissipated APE, actually occurs at the expense of 
(the exergy part of) IE. On the other hand, the APE dissipated by molecular diffusion is 
found to be dissipated into (the dead part of) IE, i.e., the same kind of IE the viscously 
dissipated KE is converted into, not into GPE r . Turbulent stirring, therefore, should 
not be viewed as introducing a new form of mechanical-to- mechanical APE /GPE r con- 
version, but simply as enhancing the existing IE/GPE r conversion rate, in addition to 
enhancing the viscous dissipation rate of KE, as well as the diffusive entropy produc- 
tion and APE dissipation rates. These results are important, for they significantly alter 
the current understanding about the nature of turbulent diapycnal mixing and its links 
with the dissipation of mechanical energy and turbulent increase of GPE r . In particular, 
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the possibility that GPE r may decrease as a result of turbulent mixing, not necessarily 
increase as is commonly thought, is to be emphasised. Moreover, the fact that the turbu- 
lent increase of GPE r is associated with an enhanced IE/GPE r conversion physically 
implies that compressible effects must be considerably larger than previously thought, 
raising fundamental questions about the possible limitations of the widely used incom- 
pressible assumption in the modelling of fluid flows at low Mach numbers, which further 
work should aim at elucidating. Finally, the present results also have implications for 
the way one should quantify the efficiency of mixing in turbulent stratified fluids, with 
new definitions for the mixing efficiency (or more accurately dissipations ratio) ^mixing 
and flux Richardson number Rf being proposed in the introduction, where they are also 
compared with existing definitions. 

A significant achievement of the extended APE framework is to allow for a more rig- 
orous and general re-derivation of MW98's result (our Eq. (3.17)) that is also valid for a 
non-Boussinesq ocean, and which results in the appearance of a nonlinearity parameter 
£ ^ 1, with MW98's results being recovered for £ = 1. The main new result here is the 
finding that the work rate done by the surface buoyancy fluxes G(APE) should be numer- 
ically comparable to W r _f orcing . This is important, because while W r j orcing is currently 
widely agreed to be large 0(0.4 TW), G(APE) has been widely thought to be negligible 
on the basis of MW98's argument that this is required by Sandstrom's "theorem". The 
result therefore demonstrates that G(APE) is as important as the mechanical forcing in 
driving and stirring the oceans, in agreement with Oort & al. (1994)'s previous conclu- 
sions. The two main consequences are that: 1) there is no reason to reject the idea that 
the oceans are a heat engine, and that the North-South heat transport is mostly the result 
of the buoyancy forcing, as was usually thought to be the case prior to MW98's study, 
e.g. Colin de Verdiere (1993); 2) the overall value of ^mixing in the oceans is likely to be 
significantly larger than the value of ^mixing — 0.2 currently used, as seems to be required 
by the large magnitude of G(APE), given that buoyancy-driven turbulent mixing has a 
significantly higher ^mixing than mechanically-driven turbulent mixing in general. Note 
that increasing ^mixing decreases the requirement on the mechanical sources of stirring, 
thereby providing a natural way to remove the apparent paradoxical result of a shortfall 
in the mechanical stirring energy that arises when ^mixing = 

0.2 is used in Eq. (3.17) 

for £ = 1. The extended APE framework may also be used to revisit the assumed impli- 
cations of Paparella & Young (2002) 's anti-turbulence "theorem", the main conclusion 
being that horizontal convection may in fact support elevated values of diapycnal mixing, 
in contrast to what is usually assumed. It seems important to stress, however, that even 
if one could prove that surface buoyancy fluxes could sustain on their own diapycnal 
mixing rates as large as observed, it would not prove that the actual buoyancy- driven 
component of the ocean circulation is not mechanically-controlled. Indeed, it is essential 
to recognise that both G(KE) and G(APE) do not depend only on the external forcing 
parameters, but also strongly depend on the actual ocean circulation and stratification, 
which necessarily implies a mechanical control of the buoyancy-driven circulation and 
vice versa. In other words, even though the present study disagrees with MW98's con- 
tention that surface buoyancy fluxes do not work or stir the oceans, there should be 
little doubt that even if the MOC is expected to be primarily buoyancy-driven, it is also 
mechanically controlled in ways that are not fully understood. In this respect, therefore, 
MW98 deserves credit for challenging the idea that the behaviour of the buoyancy- driven 
circulation can be understood independently of its link with the mechanical forcing. 

An important challenge ahead will be to extend the present framework to deal with 
an equation of state also depending on salinity, which gives rise to the possibility of 
double-diffusion effects, and storing energy in chemical form. Such extensions are needed 
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to better connect the present results to many laboratory experiments based on the use of 
compositionally-stratificd fluids, although the present evidence is that the use of a linear 
equation of state is probably accurate enough to describe the latter. To that end, many 
technical difficulties need to be overcome. Indeed, salinity complicates the definition of 
Lorenz (1955)'s reference state to such an extent that it is not even clear that such a 
state can be uniquely defined, e.g. Huang (2005), as may be the case in presence of 
humidity in the atmosphere, e.g. Taillcux & Grandpeix (2004). A potentially important 
generalisation would also be to further decompose the internal energy in order to isolate 
the available acoustic energy considered by Bannon (2004), which in the present paper is 
included as part of our definition of APE. Finally, further work is required to understand 
how to fix the value of £ and ^mixing in Eq. (3.17) in the actual oceans. 
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Appendix A. Energetics of Incompressible Navier-Stokes Equations 

A.l. Boussinesq equations with equation of state nonlinear in temperature 

The purpose of this appendix is to document the energetics of the Boussinesq system of 
equations that form the basis for most inferences about stratified turbulence for fluid flows 
at low Mach numbers, and which is commonly used in the theoretical and numerical study 
of turbulence, e.g. Winters & al. (1995); Caulfield & Peltier (2000); Staquet (2000); 
Peltier & Caulfield (2003). In order to go beyond the usual case of a linear equation of 
state, a slight generalisation is introduced by allowing the thermal expansion coefficient 
to vary with temperature. The resulting set of equations is therefore as follows: 

/)v 1 . V P=-^z + ,V 2 v (Al) 



Dt 



Po 



V- v 

DT 
~Dt 



Po 







= kV 2 T 



P(T) = Po 



1 



f 



a(T')dT' 



(A 2) 
(A3) 

(A 4) 



where v = (u, v, w) is the three-dimensional velocity field, P is the pressure, p the density, 
T the temperature, v ~ p/ p the kinematic viscosity, p the (dynamic) viscosity, k the 
molecular diffusivity, g the acceleration of gravity, and p a reference density. The classical 
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Boussinesq model, called L-Boussinesq model in this paper, is simply recovered by taking 
a to be a constant in Eq. (A 4). In that case, Eqs. (A3) and (A 4) may be combined to 
obtain the following diffusive model for density: 

^ = «V 2 P (A 5) 

as assumed in many numerical studies of turbulence, e.g. Winters & al. (1995); Caulfield 
& Peltier (2000); Staquet (2000); Peltier & Caulfield (2003). When the temperature 
dependence of a is retained, the resulting model is called here the NL-Boussinesq model. 

A. 2. Standard energetics 

Evolution equations for the KE and GPE are obtained by the standard procedure, e.g. 
Batchclor (1967); Landau & Lifshitz (1987), assuming that the system is forced me- 
chanically by an external stress r, and thermodynamically by external heat fluxes, both 
assumed to act at the surface boundary located at z = 0. The first equation is premulti- 
plied by po v an d volume-integrated. After re-organisation, the equation becomes: 

d{KE) 



dt 



= i f Po^rdV = f r-u s dS- f P gwdV+ f Po edV (A 6) 

dt Jy 2 Jgy Jy Jy 



G(KE) W D(KE) 

where W is the so-called density flux, D(KE) is the viscous dissipation rate of kinetic 
energy, and G(KE) is the rate of work done by the external stress. The time evolution 
of the total gravitational potential energy of the fluid, i.e., the volume integral of pgz, is: 

d(GPE) 



dt 



= 4 / pgzdV= [ pgwdV- [ p gzanV 2 TdV, (A 7) 

dt Jy Jy Jy 



W 



where B is the Boussinesq approximation of the work of expansion/contraction. In the 
present case, it is possible to derive an explicit analytical formula for B: 

f f dT f dcx 

B = - gzp a K V 2 TdV= / K g Po a—dV+ / K p gz—{T) || VT|| 2 dV 



V JV u ^ JV 



q / ruj" lr ri < 

oz Jy dT 



= Kg[(p) bot tom-(p)to P }+ np a gz-^(T)\\VT\\ 2 dV, (A 8) 
. ' Jv dl 

by using integration by parts, and using the fact that the surface term vanishes because 
the surface is by assumption located at z = 0, where (p)bottom and (p)top denote the 
surface-integral of the bottom and top value of density. For a linear equation of state, 
B = Bl will in general be small, because of the smallncss of the molecular diffusivity k, 
and finite top-bottom density difference. When a increases with temperature, however, B 
may become significantly larger than Bl in turbulence strong enough as to make ||VT|| 2 
large enough for the second term in Eq. (A 8) to overcome Bl, pointing out the possibly 
critical role of nonlinearity of the equation of state in strongly turbulent fluids. 

A. 3. Lorenz (1955) 's available energetics 

We now seek evolution equations for the available and un-available parts of the gravita- 
tional potential energy, as previously done by Winters & al. (1995) in the case of the 
L-Boussinesq equations. By definition, the expression for the GPE r is: 

GPE r = ( p r gz r dV, (A 9) 

Jv 
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where z r = z r (x,t) and p r = p r (z r ,t) are the vertical position and density of the parcels 
in the Lorenz (1955) 's reference state. In Boussinesq models, fluid parcels are assumed 
to conserve their in-situ temperature in the reference state, so that T r (z r ,t) = T(x,t). 
Taking the time derivative of Eq. (A 9) thus yields: 

d(GPE r ) f D Pr „, f Dz r 
= gz r 



dt 



J 9Zr^dV + gPr^dV = ~ f gz r p aK\7 2 TdV 



=0 



= - / gz r p aK,VT ■ n dS + / np gVT ■ V(az r )dV 

JdV JV 

= -W rJorcmg + J Kp g\\Vz r \\ 2 ^^^dV, (A 10) 

. 1 1 — / 

W r , mixing 

where W r j orc i ng is the rate of change of GPE r due to the external surface heating/cooling. 
For a Boussinesq fluid, this term is identical to the APE production rate G(APE), as 
shown in the following, i.e., W r j orc i ng = G(APE). The above formula was obtained by 
using the following intermediate results: 

V[T(x)] = V[T r (z r (x))] = ^Vz r , (All) 
V[a(7> r ] = V[a(T r (z r ))z r ] = ^^Vz r , (A 12) 

OZ r 

as well as the important result that the integral involving the term Dz r /Dt vanishes 
identically established by Winters & al. (1995). In this paper, the result was established 
by using an explicit formula for the reference stratification. An alternative way to recover 
such a result is achieved by noting that the velocity v r = (Dx r /Dt, Dz r /Dt) of the fluid 
parcels in the reference state must satisfy the continuity equation: 

V r • v r = 0, (A 13) 

where V r is the divergence operator in the reference space state (x r , z r ), from which it 
follows that the surface-integral of W r ^ m i X i ng = Dz r /Dt along each constant z r level must 
vanish, which implies Winters & al. (1995) s result. The equation for APE = AGPE = 
GPE - GPE r becomes: 

d(APE) _ d(GPE) _ d{GPE r ) 
dt dt dt 

= W- {Waiving - B) + W rJorcmg = G(APE) + W - D(APE), (A 14) 

where 

D{APE) = W r , mixing - B 



f n?9T r f da 

= J pogaK\\Vz r \\ —dV - B L + J Kp g{z r - z) — \\\/z r \ 



fdT r x 1 



dv > (A15) 
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by using the results that: 



/ v A»«ilVT|^-/«^|V Jp p(g)V. (A16) 
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f „_ .. 9 d(az r )dT r _ r f .._ .., / da dT r \ dT r ITr . , „ . 

i/^ii v ^ii -kr^ y = //o 5 Kiiv, r |i ^i +2 ,__j_^. ( A it) 

Empirically, it is usually found that D(APE) > 0, which is not readily apparent from 
the form of D(APE), and for which a rigorous mathematical proof remains to be es- 
tablished. Interestingly, while W r ^ m ixing and B appear to be both strongly modified by 
a temperature-dependent a, this is much less so for their difference D(APE), which is 
usually found empirically to be well approximated by its "linear" part Dl(APE). This 
is important, because it clearly establishes that D(APE) and W r ^ m i X i ng may be signifi- 
cantly different when the temperature dependence of a is retained, in contrast to what 
is generally admitted based on the L-Boussinesq model. This suggests that results based 
on the study of the L-Boussincsq model are likely to be more robust and accurate for the 
description of the KE/APE dynamics than for the description of GPE r . The condition 
for \D NL (APE)\ < \D L {APE)\ to be satisfied is that da/dT\dT r /dz r \\z r - z\ < 1, which 
appears to be satisfied in practice for water or seawater. Whether this is also true for 
other types of fluids still needs to be established. 



Appendix B. Energetics of compressible Navier-Stokes Equations 

B.l. Compressible Navier-Stokes Equations (CNSE) 

The purpose of this appendix is to generalise Winters & al. (1995) 's results to the fully 
compressible Navier-Stokes equations, which are written here in the following form: 

P^ + VP=- A9 z + V-S (Bl) 
^+pV-v = (B2) 



DE Q pe - V • F 



9 



Dt T P T ( B3 ^ 

J = /(E,u), (B4) 

T = T(^v) = — , P = P(E, «) = - — . (B5) 

In the present description, the three-dimensional Eulcrian velocity field v = (u, v,w), 
the specific volume v = l/p (with p the density) , and the specific entropy E are taken as 
the dependent variables, with the thermodynamic pressure P and in-situ temperature T 
being diagnostic variables as expressed by Eqs. (B4-B5), where I is the specific internal 
energy, regarded as a function of E and v. Further useful notations are: D/Dt = d/dt + 
(v • V) is the substantial derivative, e is the dissipation rate of kinetic energy, F q = 
—kxpCpVT is the diffusive heat flux, C p is the specific heat capacity at constant pressure, 
kx is the molecular diffusivity for temperature, g is the acceleration of gravity, and z is 
a normal unit vector pointing upward. Moreover, S is the deviatoric stress tensor: 



in classical tensorial notation, e.g., Landau & Lifshitz (1987), where Einstein's summa- 
tion convention for repeated indices has been adopted, and where Sij is the Kronecker 
delta. The parameters p and A are the shear and bulk (or volume) viscosity respectively. 
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B.2. Standard energetics 

The derivation of evolution equations for the standard forms of energy in the context 
of the fully compressible Navier-Stokes equations is a standard exercise, e.g. de Groot 
& Mazur (1962); Landau & Lifshitz (1987); Grimes (2004), so that only the final 
results are given. In the standard description of energetics, only the volume-integrated 
kinetic energy (KE), gravitational potential energy (GPE), and internal energy (IE) are 
considered, viz., 

r 2 



KE = [ p^-dV, GPE= f pgzdV IE = f pI(E,v)dV, (B7) 

JV 2 Jy Jy 

idard evolution equations are respectively given by: 

^ = -/ pgwdV + J P^dm + G(KE)-D(KE)-P a ^, (B8) 



w 



d(GPE) 
Jt 



pgwdV, (B9) 



w 



d(IE) 
dt 



= J v PQ d V- J dm = D(KE) + Q heaUng - Q coo Un g - B, (B 10) 



B 

where G(KE) is the rate of work done by the mechanical sources of energy on the 
fluid, QheaUng (resp. QcooUng) is the surface- integrated rate of heating (resp. cooling) 
due to the thermodynamic sources of energy, and V i is the total volume of the fluid; 
additional definitions and justifications are given further down the text. Summing Eqs. 
(B 8-B 10) yields the following evolution equation for the total energy total energy TE = 
KE + GPE + IE: 

—j^ - = G(KE) + QheaUng — QcooUng — Pa i (B 11) 

which says that the total energy of the fluid is modified: 

• by the rate of work done by the mechanical sources of energy; 

• by the rate of heating/cooling done by the thermodynamic sources of energy; 

• by the rate of work done by the atmospheric pressure P a against the volume changes 
of the fluid. 

As these derivations are quite standard, justifications for the above equations are only 
briefly outlined. Thus, the KE equation (B8) is classically obtained by multiplying the 
momentum equation by v, and integrating over the volume domain. The term W results 
from the product of v with the gravitational force vector, whereas the product v • VP = 
V-(Pv) — PV-v = V-(Pv) — (P/v)Dv/Dt yields the work of expansion/contraction minus 
the work done by the atmospheric pressure against total volume changes. The product of 
the velocity vector with the stress tensor is written as the sum G(KE) — D(KE), where 
G(KE) represents the work input due to the external stress, and D(KE) the positive 
dissipation of kinetic energy. The general expression for the mechanical energy input is: 

G{KE) = [ vS • ndS = [ t ■ vdS (B 12) 

JdV JdV 

where vS is the vector with components (Sv)j — Sijm, while Sn = t is the stress 
applied along the surface boundary enclosing the fluid. G(KE) is therefore the work 
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of the applied stress done against the fluid velocity. If one assumes no-slip boundary 
condition on all solid boundaries, then this work is different from zero only on the free 
surface. The function D{KE) is the dissipation function: 

««*>-fM% + %-&%)' +, * vf h <Bi3) 

where again the summation convention for repeated indices has been used, e.g., Landau 
& Lifshitz (1987). The equation for GPE (B9) is simply obtained by taking the time 
derivative of its definition, using the fact that D(pgzdV)/ Dt = pgw, since D(pdV)/ Dt = 
from mass conservation. The equation for IE (B 10) results from the fact that the 
differential of internal energy in the entropy/specific volume representation is given by 
dl = TdYj — Pdv. The terms Qheating and Q coo iing represents the surface-integrated net 
heating and cooling respectively going through the surface enclosing the domain. 

B.3. Available energetics 

In this paragraph, we seek to derive separate evolution equations for the available and un- 
available parts of the total potential energy PE = IE+GPE+P a V a i, as initially proposed 
by Lorenz (1955), building upon ideas going back to Margulcs (1903). Specifically, PE 
is decomposed as follows: 

PE= [ p [/(£, v) + gz] dV + P a V ol 
Jv 

= [ p[I(£,v) + gz] dV - [ p[I(Z,v r )+gz r ] dV + P a (V ol - V ol , r ) 
Jv Jv 

S v ' 

APE 

+ f p [I(S, v r ) + gz r ] dV + PaV oLr (B 14) 

Jv 

V v ' 

PE r 

where PE r is the potential energy of Lorenz (1955)'s reference state, and APE = 
PE — PE r is the available potential energy. As is well known, the reference state is the 
state minimising the total potential energy of the system in an adiabatic re-arrangement 
of the fluid parcels. From a mathematical viewpoint, Lorenz (1955) 's reference state can 
be defined in terms of a mapping taking a parcel located at (x, t) in the given state to 
its position (x r , t) in the reference state, such that the mapping preserves the specific 
entropy S and mass pdV of the parcel, viz., 

S(x, t) = S(x r , t) = Y, r (z r , t), (B 15) 

p(x, t)dV — p(x r , t)dV r — p r {z ri t)dV r , (B 16) 

where the second condition can be equivalently formulated in terms of the Jacobian 
J = <9(x r )/<9(x) of the mapping between the actual and reference state as follows: 

p(x, t) = p(x r , t)^= p r (z r , t)^. (B 17) 

Prior to deriving evolution equation for PE r and APE, it is useful to mention three 
important properties of the reference state, namely: 

(a) The density p r — p r (z r ,t) and pressure P r — P r (z r ,t) of the Lorenz (1955)'s 
background reference state are functions of z r alone (and time); 
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(b) The background density p r and pressure P r are in hydrostatic balance at all times, 
i.e., dPr/dzr = —p r g (this is a consequence of the reference state being the state min- 
imising the total potential energy in an adiabatic re-arrangement of the parcels); 

(c) The velocity v r = (Dx r /Dt, Dy r /Dt, Dz r /Dt) of the parcels in the reference state 
satisfies the usual mass conservation equation: 

Dv 

-^= Wr V r -v r , (B18) 

where V r • v r is the velocity divergence expressed coordinates system of the reference 
state, which is a consequence of the mass of the fluid parcels being conserved by the 
mapping between the actual and reference state. 

Eq. (B 18) is important, for it allows an easy demonstration of the following result: 

[ pPr^dV= f p r P r ^dV r = f P r V r -W r dVr 

Jv ut Jv r ut Jv r 

= f P r v r -n r dS r — f v r ■ VP r dV r — P a '^'° l ' r + / p r gw r dV r , 
Jdv r Jv r dt J Vr 



(B19) 



w r 

which establishes the equivalence between the work of expansion and the work against 
gravity in the reference state, where n r is a outward pointing unit vector normal to the 
boundary dV r enclosing the fluid in the reference state. In Eq. (B 19), the first equality 
stems from expressing the first integral in the reference state; the second equality uses 
Eq. (B18); the third equality results from an integration by parts; the final equality 
stems from that P r depends on z r and t only, and that it is in hydrostatic balance, and 
from using the boundary condition v r • n r = w r = dt] r /dt at the surface assumed to be 
located at z r = f] r (t). 

B.4. Evolution of the background potential energy PE r 

We seek an evolution equation for the background PE r by taking the time derivative of 
the expression in Eq. (B 14), which yields: 



d(PE r ) 
Jt 
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DE Dv r 

T r — P r — — h qw r 

Dt Dt y 



„r r, dVoi r f . T/ ^ . . D(pdV) 

+ Pa ^f + Jy [ ( ' ^ + 



Dt 



=o 



= J fT T ^dV = J pQdV + J p{^j^QdV 

= Qnet + (1 - le)D(KE) + D(APE) - G(APE), (B 20) 

where the final result was arrived at by making use of Eq. (B19), as well as of the 



definitions: 



/ P QdV= f {V • {npCpVT) + pe} dV = Q ne t + D(KE), (B21) 
Jv Jv 

Qnet — / KpC p VT ■ ndS = Qheating — Qcooling (B 22) 

Js 

G(APE)= J { ^ T Tr ^j npC p \7T-ndS (B 23) 
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Figure 10. Successive refinements of the energetics of a forced/dissipated stratified fluid. Panel 
(I): KE I APE I PE r representation. Panel (II): Decomposition of PE r into lE r + GPE r . Panel 



(III): Decomposition of IE r into a dead part IEo and exergy part IE exe 



IEr-IE () . Panel 



(IV): Decomposition of APE into AIE and AGPE, revealing the link between C(KE,APE) 
to the density flux W and work of expansion/contraction B. 



D(APE) = J k P C p VT • V (^-jr^j dV, (B 24) 

le D(KE) = J v (^Y^j pe dV (B 25) 

where n is the unit normal vector pointing outward the domain. Eq. (B 22) expresses the 
net diabatic heating Q ne t due to the surface heat fluxes as the sum of a purely positive 
Qheating and negative —Qcooiing contributions. Eq. (B 23) defines the rate of available 
potential energy produced by the surface heat fluxes. The term D(APE), as defined 
by Eq. (B24), is physically expected to represent the rate at which APE is dissipated 
by molecular diffusion, so that it is expected to be positive in general, which has been 
so far only established empirically using randomly generated temperature fields, but a 
rigorous mathematical proof is lacking. Finally, Eq. (B 25) states that a tiny fraction of 
the diabatic heating due to viscous dissipation might be recycled to produce work. If "f e 
could be proven to be positive, it could probably be included as part of the G(APE). In 
the following, it will just be neglected for simplicity. 

B.5. Evolution of Available Potential Energy (APE) 

In the previous section, we defined the total potential energy as the sum of GPE, IE, 
and the quantity P a V Q i, see Eq. (B14). As a result, using the evolution equations for 



Energetics and thermodynamics of turbulent molecular diffusive mixing 37 
GPE and IE previously derived, the evolution equation for PE is given by: 

d -m =W _ B + Q net + Pa ^L. (B 26) 
dt dt 

Now, combining this equation with the one previously derived for PE r allows us to derive 
the following equation for the available potential energy APE = PE — PE r , defined as 
the difference between the potential energy and its background value: 

d ^ t E) ~ W - B + Pa ^ d T +Qnet+D(KE)- [D(KE) + D(APE) + Q net - G(APE) 



C(KE,APE) 



= C{KE, APE) + G(APE) - D(APE) , (B 27) 

where the final expression neglects the small term j e D(KE). The corresponding energy 
flowchart for the KE / APE / PE r system is very simple, and is illustrated in Fig. 10 (Panel 
I). This diagram shows that mechanical energy enters the fluid via the KE reservoir, 
and that thermal energy enters it via the PE r reservoir. There are two dissipation routes 
associated with the viscous dissipation of KE and the diffusive dissipation of APE. Only 
a certain part G(APE) of the thermodynamic energy input can be converted into APE 
and hence into KE, which is processed via the PE r reservoir. The two-headed arrow 
indicates the reversible conversion between KE and APE. 



B.6. Splitting of PE r into GPE r and IE r 

Although the KE / APE / PE r system offers a simple picture of the energetics of a (turbu- 
lent or not) stratified fluid, it is useful to further decompose the background PE r reservoir 
into its GPE and IE component, in order to establish the link with the existing liter- 
ature about turbulent mixing, as well as with Munk & Wunsch (1998) 's theory. The 
particular question to be addressed is to understand how much of D(KE) and D(APE) 
are actually spread over GPE r and IE r . Likewise, to what extent do Qnet and G(APE) 
affect GPE r compared to IE r , where we have the following definitions: 

IE r = I p(x,t)I{i:,v r )dV = / p r (z r ,t)I(£,v r )dV r (B28) 
Jv Jv r 

GPE r = / p(x,t)gz r (xL,t)dV = / p r (z r , t)gz r dV r , (B 29) 

Jv Jv r 

by expressing the integrals in the coordinate system associated either with the actual 

state or reference state. By definition, 

d(GPE r ) f Dz r „ T f D(pdV) f „ T TTr 

S v ' 

=0 

so that the evolution equation for IE r + P a V ^ r = PE r — GPE r is simply: 

d{IE r +P a Vol,r) 



dt 



Q net + D(KE) + D(APE) - G(APE) - W r (B 31) 



In order to make progress, we need to relate W r to the different sources and sinks affecting 
PE r , as identified in Fig. (10). To that end, we use the fact that W r is related to the 
work of expansion in the reference state, as shown by Eq. (B 19), and regard v = P) 
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as a function of entropy and pressure, for which the total differential is given by: 

dv = TdY>- -^dP (B32) 
p 2 cj 

where T = aT/(pC p ) is the so-called adiabatic lapse rate, e.g. Feistel (2003), and c 2 s = 
(dP/dp)z is the squared sound of speed. As a result, the expression for W r becomes: 



a r T r Q 1 DP r 



p r C pr T p2 c 2 r Df 



dVr (B33) 



where P' r = P r — P a is the pressure corrected by the atmospheric pressure, by noting 
that we have: 

dV i^ r f Dv r f Dv r 

— = J v ^r pdv = J Vr -w prdVr - (B 34) 

In order to simplify Eq. (B33), let us recall that mass conservation can be rewritten in 
hydrostatic pressure coordinates as follows: 

Q BP 

V,.u r + — -^=0 (B35) 

e.g. Haltiner & Williams (1980); de Szoeke & Samelson (2002). As a result, it follows that 
integrating Eq. (B 35) from the surface where P a = est, and hence where DP r /Dt = 0, 
to an arbitrary level indicates that the surface integral of DP r /Dt must vanish along any 
z r = est surfaces. As a consequence, the term depending on DP r /Dt in Eq. (B33) must 
vanish. For an alternative derivation of this result, see Pauluis (2007). The remaining 
term can be written as follows: 

W r = I ^f^{V-(*pC P VT) + pe} dV 

Jy Pr^pr 1 

= i ^(i + ^)v.( K pC p vr) dV+j ^% P edv 

Jy p r \j V r \ 1 J Jy Pr^pr 1 

I ^V.(^VT) dV+j ^(Tl-L\ V . {Kp CpVT) dV+f ^% P edV 

Jy Pr^pr Jv Pr^pr \ 1 / JV Pr^pr 1 

s v ' s v ' s v ' 

W T , mixing -W r j orcing T rapc D(APE) T rikc D(KE) 

= W^mixing — W r j orc i ng + T r ^ ape D{APE) + T r ^ e D(KE) (B36) 
where we defined: 

W„ mg = - f npCpVT ■ V (^-) dV (B 37) 

JV \Pr^pr J 

r rv P' 

W rJorcing = - / -J—^npC p VT -ndS (B 38) 

is Pr^pr 

T r>ape D{APE) = f (^X) V • (npCpVT) dV (B 39) 

Jy p r (spr \ 1 J 
I" cv P T 

T rM D{KE)= ^-^^pedV (B40) 

JV Pr^pr 1 

Eq. (B 36) shows that the variations of GPE r are affected by: 

• Turbulent mixing, associated with W r fixing- This expression is similar to the one 
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previously derived for the L-Boussinesq model. The classical Boussinesq expression can 
be recovered from using the approximation T w T r , taking a r , p r , and C pr as constant, 
and using the approximation P r w= —p gz r , which yields: 

f - dT r 

W r , mixing « J K,p g\\Vz r \\ a-^-dV; 

• The surface forcing, associated with W r j orc i ng . Likewise, the L-Boussinesq expres- 
sion can be recovered by making the same approximation, yielding: 

w ~ I agZr n A Q 

V* r,j 'orcing ~ / ^ tysurfU'D 
JV ^p 

Note that in the L-Boussinesq approximation, we have: 

W rJorcmg « G(ALL) 

which is not generally true in the fully compressible Navier-Stokes equation; 

• The contribution from the viscous and diffusive dissipation of KE and APE respec- 
tively associated with D(KE) and D(APE). Note that the coefficient T r , ap e and T r .ke 
are very small for a nearly incompressible fluid such as seawater. For instance, typical 
values are: a = 10~ 4 K-\ P = 4.10 7 Pa, p = 10 3 kg. m - 3 , and C p = 4.10 3 J.kg-^K" 1 , 
which yields: 



V 10 3 x 4.10 3 J y > 



From this, it follows that at leading order, the direct effects of D(APE) and D(KE) on 
GPE r can be safely neglected compared to the other two effects, so that: 

d(GPE ) 

^ Wr ~ W^r,mixing Wr, forcing (-^^-O 

The resulting modifications to the energy flowchart are then displayed in Fig. 10 (Panel 
II). At leading order, the effects of the forcing and mixing on GPE r appear as conversion 
terms with IE r . 

B.7. Further partitioning of internal energy into a "dead" and "exergy" component 
As seen previously, the L-Boussinesq model is such that: 

D(APE) w W r , mixing , G{APE)nW rJorcing , (B42) 

which may give the impression, based on Fig. 10 (Panel II) that the APE dissipated by 
D(APE) is actually converted into GPE ri while G(APE) may also appear as originating 
from GPE r . The purpose of the following is to show that this is actually not the case. 
To that end, we introduce an equivalent isothermal state having exactly the same energy 
as Lorenz (1955)'s reference state, i.e., that is defined by: 

IE r + GPE r + P a V ol , r = I E + GPE + P a V olfi . (B 43) 

^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^ 

PE r PEo 

Because both Lorenz 's reference state and the equivalent thermodynamic equilibrium 
state are in hydrostatic balance at all time, PE r and PE are just the total enthalpies of 
the two states. This makes it possible to define each parcel by their horizontal coordinates 
(x, y) and hydrostatic pressure P, and to assume that the dead state can be obtained 
from Lorenz's reference state by an isobaric process, so that (xo,yo,Po) = {x ri y ri P r ), 
which in turn implies (dxo/dt,dyo/dt,d,Po/dt) = (dx r / 'dt, dy r / 'dt, dP r / 'dt) . 
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Prior to looking at the evolution of the dead state, let us establish that if the pressure 
P is in hydrostatic balance at all times, then we have the following result: 

J J^dV = J u-V h PdV, (B44) 

where u is the horizontal part of the 3D velocity field, and Vh the horizontal nabla 
operator. The proof is: 

= j t \PaVoi + M tot gH b + pgzdV^ - pj¥j± - 1. j pgzdV = 
where V a i and M tot are the total volume and mass of the fluid, whose expressions are: 

Vol = I {n{x,y,t) + H b )dxdy, gM tot = / (P b (x,y,t) - P a )dxdy, 
Js Js 

where z = r/(x,y,t) is the equation for the free surface, P b (x,y,t) is the bottom pres- 
sure, H b is the total depth of the basin, and where the expression between brackets was 
obtained by using the following result: 

f PdV = [ [Pz] v _ H dxdy+ [ pgzdV = [ [P aV + P b H b ]dxdy + [ pgzdV 
Jv Js Jv Js Jv 

= P a [ (r] + H)dxdy + H [ (P b - P a )dxdy + [ pgzdV. 
Js Js J 

S v ' S v ' 

P a V ol M tot gH 

The important consequence of Eq. (B 44) is that the volume integral of DP/ Dt identically 
vanishes when P is independent of the horizontal coordinates, as is the case for P r and 
Po- Now, using the expression for the enthalpy I + P/p: 

d(I + P/p) = C p dT + (v- T^j dP = C p dT + v(l- uT) dP, 

we can derive the following equation for PEq, 

d(PE ) f ( DT a DP \ dT f 

-J Pa y°p0^f + U - a o T o) -jyj-J dV = —J p C p0 dVo, (B45) 



dt 

which naturally provides the following equation for T : 

dT _ D(KE) + D(APE) + Q- G(APE) 
dt J Vo p CpodVo 

We can now derive an evolution equation for GPE , using the relation: 

d(GPEo) 



(B46) 



= J Po^PodVo (B47) 



dt j Vo 

where Pq = Pq — P a . Now, expressing dv — vadT + wydP, where 7 is the isothermal 
expansion coefficient, we arrive at the following expression: 



d{GPE ) 
dt 



PoPa 



DT Q DP 
vo—pr: v Jo 



Dt ' Dt 



dV = ^ J aoP^dVoPoCpodVo (B48) 
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noting again that the term proportional to DPo/Dt must vanish from the arguments 
developed above, so that we simply have: 

d{GPE ) 



dt 

where 



= T \D{KE) + D(APE) + Q net - G{APE)\ (B 49) 

Y Ivo Po a odVo ^ 
Iv PoCpodVo 

As a result, it follows that: 

djlEp + P a V olfl ) = d{PE - GPE Q ) 
dt dt 

= (1 - T ) [D{KE) + D(APE) + Q net - G(APE)\ (B 51) 

Let us now define the exergy part of the IE r + P a Voi,r as 

IE exer gy = I E r — I Eq + P a (V ^ r — V lfl) (B 52) 

The equation is: 



d{l E eXer gy) 



~W r + T 



D{KE) + D(APE) + Q net - G(APE) 



dt 

= (T - T rM )D{KE) + (T - T r%ape )D(APE) + T [Q net - G(APE) 

^^r,mixing ~t~ ^^r,forcing- (B 53) 

Again, the neglect of the terms proportional to aP/(pC p ) yields the following simplifi- 
cation: 

d(IE + P a V olfi ) 



D(APE) + D(KE) + Q net - G(APE), (B 54) 



dt 

d(IE exergy ) ^ _w i'r^^'i 

"V r.rnixing *v r.f orcing • yr> Do ) 



dt 

The corresponding energy flowchart is this time illustrated in Fig. 10 (Panel III). This 
figure shows that when IE r is decomposed into its dead and exergy part, a decoupling 
between the KE/APE/IE and IE r — IE /GPE r reservoirs appears at leading order. 
Note, however, that the rates between the reservoirs remain coupled, owing to the cor- 
relation between D(APE) and W r ^ m i X i n g, as well as between G(APE) and W r j orc i ng 
discussed in this paper, and which is a central topic of turbulent mixing theory. 

B.8. Separate evolution of APE into GPE and IE components 

We conclude the evolution equation of energetics by further splitting the APE reservoir 
into its GPE and IE components. Using the previous relations, one easily shows that: 

d(AGPE) 



dt 

d(AIE) 



= W - W r W W - W r , mixing + W r J orc i n g (B 56) 



W - B + G(APE) - D(APE) - [W - W„ mg + W rJorcmg } 

~ —B + W rtTnixing - D(APE) + G(APE) - W r j orcmg . (B 57) 

For seawater, it is generally found that AIE accounts for around 10% of the total 
APE, so that to a good approximation APE w AGPE, which is implicit in the Boussi- 
nesq approximation. Equating d(AGPE)/dt with d(APE)/dt amounts to requiring that 
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d(AIE)/dt w 0. By imposing that the forcing and mixing terms vanish separately, one 
obtains: 

D(APE) w W r , mixmg - B, (B 58) 

G(APE) « W r 

, forcing •> 

(B59) 

which are equivalent to those of the L-Boussinesq and NL-Boussinesq model. The cor- 
responding energy flowchart is depicted in Fig. fO (Panel IV). The key feature of this 
figure is to reveal that the conversion rates between AGPE and AIE are identical to 
those taking place between IE r — IEo and GPE r , which appears to be where the coupling 
between stirring and mixing fundamentally occurs. 
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